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Bulge n and B/T in High Mass Galaxies: Constraints on the 
Origin of Bulges in Hierarchical Models 

Tim Weinzirl^, Shardha Jogee^, Sadegh Khochfar^'^, Andreas Burkerf*, John Kormendy^ 

ABSTRACT 

We use the bulge Sersic index n and bulge-to-total ratio {B/T) to explore 
the fundamental question of how bulges form. We perform 2D bulge-disk-bar 
decomposition on H-hand images of 143 bright, high mass {M^ > 1.0 x IO^^Mq) 
low-to-moderately inclined {i < 70°) spirals. Our results are: (1) Our i7-band 
bar fraction (~ 58%) is consistent with that from ellipse fits. (2) 70% of the 



^ . stellar mass is in disks, 10% in bars, and 20% in bulges. (3) A large fraction 

(^ 69%) of bright spirals have B/T < 0.2, and ~ 76% have low n < 2 bulges. 

These bulges exist in barred and unbarred galaxies across a wide range of Hubble 
^ : types. (4) About 65% (68%) of bright spirals with n < 2 (B/T < 0.2) bulges 

host bars, suggesting a possible link between bars and bulges. (5) We compare 
fr^ . the results with predictions from a set of ACDM models. In the models, a high 

^ , mass spiral can have a bulge with a present-day low B/T < 0.2 only if it did not 

"Tj" ! undergo a major merger since z < 2. The predicted fraction (~ 1.6%) of high 

Q ! mass spirals, which have undergone a major merger since z < A and host a bulge 

h^ ! with a present-day low B/T < 0.2, is a factor of over thirty smaller than the 

O I observed fraction (~ 66%) of high mass spirals with B/T < 0.2. Thus, contrary 

Q ! to common perception, bulges built via major mergers since z < 4 seriously fail 

to account for the bulges present in ~ 66% of high mass spirals. Most of these 
k> ! present-day low B/T < 0.2 bulges are likely to have been built by a combination 

H I of minor mergers and/or secular processes since z < 4. 

Subject headings: galaxies: bulges — galaxies: evolution — galaxies: forma- 
tion — galaxies: fundamental parameters — galaxies: interactions — galaxies: 
structure 
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Introduction 



The formation of galaxies is a classic problem in astrophysics. Contemporary galaxy 
formation models combine the well-established A Cold Dark Matter (ACDM) cosmology, 
which describes behavior of dark matter on very large scales, with baryonic physics to model 
galaxy formation. In the early Universe, pockets of dark matter decoupled from the Hubble 
flow, collapsed into virialized halos, and then clustered hierarchically into larger structures. 
Meanwhile, gas aggregated in the interiors of the halos to form rotating disks, which are 
the building blocks of galaxies (Steinmetz & Navarro 2002; Cole et al. 2000). Such disks 
are typically destroyed during major mergers of galaxies with mass ratio M^jM^ > 1/4 (e.g. 
Steinmetz & Navarro 2002; Naab & Burkert 2003; Burkert & Naab 2004; Mihos & Hernquist 
1996). When the mass ratio is close to unity, the remnant is a spheroid with properties close 
to that of a classical bulge, namely a steep de Vaucouleurs r^/^ surface brightness profile 
and a high ratio of ordered-to-random motion {v/a). We shall return to this point in § [51 
Within this hierarchical framework, the disk of spiral galaxies forms when gas of higher 
specific angular momentum subsequently accretes around the bulge (Steinmetz & Navarro 
2002; Burkert & Naab 2004). 

ACDM-based simulations of galaxy formation face several challenges. One issue is the 
angular momentum problem; simulated galaxy disks have smaller scalelengths and, therefore, 
less specific angular momentum than their counterparts in nature (Navarro & Steinmetz 2000; 
Burkert & D'Onghia 2004; D'Onghia et al. 2006). A second issue is the problem of bulgeless 
or low bulge-to-total mass ratio {B/T) spirals. Within the ACDM paradigm, galaxies that 
had a past major merger at a time when its mass was a fairly large fraction of its present- 
day mass are expected to have a significant bulge with large B/T and high Sersic index. 
Depending on the merger history and hence the fraction of spiral galaxies that fulfill this 
criterion (see § 15.81) we can end up with a small or large fraction of present-day galaxies with 
low B/T. 

There is rising evidence that low B/T and bulgeless galaxies are quite common in the 
local Universe, especially in low mass or late-type galaxies. Late-type Sd galaxies often 
harbor no bulge (Boker et al. 2002). Kautsch et al. (2006) and Barazza, Jogee & Marinova 
(2007, 2008) also find from the analysis of several thousand late-type SDSS galaxies that 
15-20% of such disk galaxies out to z ~ 0.03 appear bulgeless. Of the 19 local galaxies 
(D < 8 Mpc) with circular velocity Vc > 150 km s~^, 11 (58%) have pseudobulges instead 
of merger-built classical bulges (Kormendy & Fisher 2008). 

Theoretical work by Koda et al. (2007) conclude the survival of disk-dominated systems 
in a ACDM universe is compatible with observational constraints provided classical bulges 
form only in mergers where M1/M2 > 0.3 and the primary halo has virial velocity V^r > 
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55 km s ^. 

Evidence also suggests that bulges with low B/T and low Sersic index n may be common 
even in high mass and/or early-type spirals. Balcells et al. (2003) report that early-type disk 
galaxies tend to have n < 3 and often from 1 to 2. Laurikainen et al. (2007) find barred 
and unbarred early-type disk galaxies to have mean B/T between 0.25 and 0.35, while later 
Hubble types have B/T < 0.2; they also find mean bulge Sersic indices to be ~ 2.5 or 
less across the Hubble sequence. Graham & Worley (2008) report low B/D ratios across 
the Hubble sequence based on bulge-disk decomposition of K-hand images of local spiral 
galaxies. They suggest that these low values are problematic for ACDM simulations, but no 
quantitative assessment of the extent of the problems is presented. 

These emerging statistics on the fraction of bulgeless {B/T ~ 0) galaxies, and galaxies 
with low B/T and low n bulges provide important first constraints. More work is needed to 
fully explore the the distribution of bulge properties in both high and low mass galaxies. In 
particular, we need to explore how the observed distributions of bulge B/T and n compare 
with the predictions from ACDM-based simulations of galaxy evolution. To the best of 
our knowledge, few such quantitative comparisons have been attempted, so that it remains 
unclear how serious the problem of low B/T galaxies is. This study is an attempt to derive 
robust observational constraints on bulge properties in high mass spirals and to attempt 
such a comparison with models. 

Completely resolving the issue of low B/T systems will require understanding the differ- 
ent types of bulges and their formation pathways. Bulges are commonly divided in several 
groups: classical bulges, boxy /peanut bulges, and 'pseudobulges' or disky bulges. Classi- 
cal bulges are believed to be built by major mergers (M1/M2 > 1/4) and the associated 
violent relaxation of stars. They are associated with modest-to-high bulge Sersic indices, 
in the range 2 < n < 6 (Hopkins et al. 2008; Springel et al. 2005; Robertson et al. 2006; 
§ 15.81) . Boxy/peanut bulges are believed to be the result of vertical resonances and buckling 
instabilities in bars, which are viewed at high inclination (Combes & Sanders 1981; Combes 
et al. 1990; Pfenniger & Norman 1990; Bureau & Athanassoula 2005; Athanassoula 2005; 
Martinez- Valpuesta et al. 2006). Pseudobulges or disky bulges are believed to form as a result 
of gas inflow into the central kiloparsec and subsequent star formation building a compact 
disky, high v/a stellar component (Kormendy 1993; Jogee 1999; Kormendy & Kennicutt 
2004, hereafter KK04; Jogee, Scoville, & Kenney 2005; Athanassoula 2005; Kormendy & 
Fisher 2005). Pseudobulges tend to have a bulge n < 2.5 (Kormendy & Fisher 2005; Fisher 
& Drory 2008). 

One possibility for the formation of disky bulges or pseudobulges is the idea of secular 
evolution (Kormendy 1993; KK04; Jogee, Scoville, & Kenney 2005), where a stellar bar or 
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globally oval structure in a non-interacting galaxy drives the gas inflow into the inner kpc via 
shocks and gravitational torque. Another idea for building disky bulges is that the gas inflow 
into the inner kiloparsec is driven by external non-secular processes, such as tidal interaction 
and minor mergers. The gas inflow in such cases can be caused by a tidally induced non- 
axisymmetric feature, such as a bar (e.g., Quinn et al. 1993; Hernquist & Mihos 1995), and 
by tidal torques from the companion. The subsequent central star formation can still form 
a compact high v/a stellar component, aka a pseudobulge. 

Throughout this paper, we avoid making any a priori assumptions about the origin of 
different types of bulges by simply referring to them according to their bulge Sersic index 
n or bulge-to-total mass ratio (B/T). We consider bulges of high (n > 4), intermediate 
(2 < n < 4) and low {n < 2) index, as well as those of low or high B/T. 

The structural properties of galaxy components, such as bulges, disks, and bars can be 
derived through the decomposition of the 2D light distribution, taking into account the PSF. 
Many early studies have performed only two component 2D bulge-disk decomposition (e.g., 
Allen et al. 2006; Byun & Freeman 1995; de Jong 1996; Simard 1998; Wadadekar et al. 1999), 
ignoring the contribution of the bar, even in strongly barred galaxies. However, recent work 
has shown that it is important to include the bar in 2D decomposition of barred galaxies, 
else the B/T ratio can be artiflcially inflated, and bulge properties skewed (e.g., Laurikainen 
et al. 2005, 2007). Furthermore, since most (> 60%) bright spiral galaxies are barred in the 
NIR (Eskridge et al. 2000; Laurikainen et al. 2004; Marinova & Jogee 2007, hereafter MJ07; 
Menendez-Delmestre et al. 2007), the inclusion of the bar is quite important. This has led 
to several recent studies, where 2D bulge-disk-bar decomposition are being performed (e.g. 
Laurikainen et al. 2007; Reese et al. 2007; Gadotti & Kauffmann 2007). 

Another advantage of bulge-disk-bar decomposition over bulge-disk decomposition is 
that the former allows us to constrain the properties of the bar itself. Bars provide the most 
important internal mechanism for redistributing angular momentum in baryonic and dark 
matter components (e.g. Weinberg 1985; Debattista & Sellwood 1998, 2000; Athanassoula 
2002; Berentzen, Shlosman, & Jogee 2006). They efficiently drive gas inflows into the central 
kpc, feed central starbursts (Elmegreen 1994; Knapen et al. 1995; Hunt & Malakan 1999; 
Jogee et al. 1999; Jogee, Scoville, & Kenney 2005; Jogee 2006) and lead to the formation of 
disky or pseudobulges (see above). Furthermore, the prominence of strong bars out to 2; ~ 1 
over the last 8 Gyr (Jogee et al. 2004; Sheth et al. 2008) suggest that bars have been present 
over cosmological times and can shape the dynamical and secular evolution of disks. Thus, 
quantifying bar properties, such as the fractional light and mass ratio (Bar/T), can yield 
insight into these processes. 

In this paper, we constrain the properties of bulges and bars along the Hubble sequence. 
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and compare our results with ACDM-based simulations of galaxy evolution. In § |2l we define 
our complete sample of 143 bright {Mb < —19.3) low-to-moderately inclined {i < 70°) 
spirals from the Ohio State University Bright Spiral Galaxy Survey (OSUBSGS; Eskridge et 
al. 2002), which is widely used as the local reference sample for bright spirals by numerous 
studies (e.g., Eskridge et al. 2000; Block et al. 2002; Buta et al. 2005; MJ07 ; Laurikainen 
et al. 2004, 2007). In § [3], we perform 2D bulge-disk and bulge-disk-bar decompositions of 
if -band images using GALFIT (Peng et al. 2002), and derive fractional light ratios {B/T, 
Bar/T, Disk/T), as well as Sersic indices and half light radii or scale lengths. Tests to verify 
the robustness of our decompositions are presented in § HI In § [5|, we present our results. 
Specifically, the total stellar mass present in bulges, disks, and bars is calculated (§ 15. 2p . 
In § 15. 3[ the distribution of bulge Sersic index n and B/T as a function of galaxy Hubble 
type and stellar mass is presented, and the surprising prevalence of bulges with low Sersic 
index n and low B/T is established. A comparison with other works is presented in § 15.41 
We examine how Bar/T and bar fraction (the fraction barred disks) change as a function of 
host galaxy properties in § 15.51 In § 15.81 we compare our observed distribution of bulge B/T 
and n in high mass (M^ > 1.0 x IO^'^Mq) spirals with predictions from ACDM cosmological 
semi-analytical models. § E] summarizes our results. 



2. Sample Properties 

2.1. OSUBSGS 

Our dataset is derived from the 182 if-band images from the public data release of 
the Ohio State University Bright Spiral Galaxy Survey (OSUBSGS; Eskridge et al. 2002). 
These galaxies are a subset of the RC3 catalog that have rriB < 12, Hubble types < T < 9 
(SO/a to Sm), D25 < 6'.5, and -80° < 5 < +50°. Imaging of OSUBSGS galaxies spans 
optical and near infrared (NIR) wavelengths with BVRJHK images available for most 
galaxies. OSUBSGS images were acquired on a wide range of telescopes with apertures 
ranging from 0.9-2. 4m. The JHK data were acquired with a variety of telescopes and 
detectors, but mainly with the 1.8 m Perkins reflector at Lowell Observatory and the CTIO 
1.5 m telescope with the OSIRIS detector, having 18.5 micron pixels (Depoy et al. 1993). 
Pixel scale is dependent on the telescope and for these observations ranged between 1 — 
1.50"/pix. Exposure times were heterogeneous, but the total observing time per object was 
typically between 10-15 minutes in H . The resulting limiting if-band surface brightness 
is ~ 20 mag arcsec"^. The typical limiting surface brightnesses of the images ~ 26 mag 
arcsec"^ in 5-band and ~ 20 mag arcsec"^ in if -band (Eskridge et al. 2002). The seeing 
depends on observing time and location. We find the if -band images have seeing of ~ 3". 
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We choose to use the NIR images rather than optical ones for several reasons. Firstly, 
NIR images are better tracers of the stellar mass than optical images, and the mass-to-light 
ratio is less affected by age or dust gradients. Secondly, obscuration by dust and SF are 
minimized in the NIR, compared to the optical. As the i^-band images are of poor quality, 
we settle on using the if -band images. 

The OSUBSGS is widely used as the local reference sample for bright spirals by numer- 
ous studies (e.g. Eskridge et al. 2000; Block et al. 2002; Buta et al. 2005; MJ07; Laurikainen 
et al. 2004, 2007). Thus, there are numerous complementary results that we can use or com- 
pare with. In particular, MJ07 have identified bars in this sample using quantitative criteria 
based on ellipse fitting, and characterized their sizes, position angles, and elhpticities. 

OSUBSGS is a magnitude-limited survey with objects whose distances range up to 
~ 60 Mpc. Faint galaxies are inevitably missed at larger distances, resulting in the absolute 
magnitude distribution in Figure [H We compare the 5-band LF of this sample with a 
Schechter (Schechter 1976) LF (SLF) with $* = 5.488 x lO'^ Mpc'^, a = -1.07, and 
M^ = —20.5 (Efstathiou, Ellis & Peterson 1988) in Figure [21 The volume used to determine 
the number density in each magnitude bin is 

V^a. = yC?L.(M), (1) 

where 

dr^aAM) = 10l+0-2(™-*^) (2) 

is the maximum distance out to which a galaxy of absolute magnitude M can be observed 
given the cutoff magnitude rric. If the SLF is representative of the true LF, then Figure [2] 
suggests that the OSUBSGS sample is seriously incomplete at Mb > —19.3, while at the 
brighter end (-19.3 to -23) the shape of its LF matches fairly well the SLF. We thus conclude 
that the sample is reasonably complete for bright {Mb < —19.3 or Lb > 0.33 L*) galaxies. 

We exclude highly inclined {i > 70°) galaxies for which structural decomposition does 
not yield accurate results. Thus, our final sample SI consists of 143 bright {Mb < —19.3) 
low-to-moderately inclined {i < 70°) spirals with Hubble types mainly in the range SO/a to 
Sc (Figure [1]). Of the 126 for which we could derive stellar masses (see § I2.2[ most have stellar 
masses M^ > 1.0 x lO^^M© (Figure [3]). Table [1] summarizes the morphologies, luminosities, 
and stellar masses of the sample. Note that there are few galaxies of late Hubble types (Scd 
or later) and we do not draw any conclusions on such systems from our study. In a future 
paper, we will tackle galaxies of lower mass and later Hubble types. 
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2.2. Stellar Masses 

We derive global stellar masses for most of the OSUBSGS sample galaxies using the 
relation between stellar mass and rest-frame B — V color from Bell et al. (2003). Using 
population synthesis models, the latter study calculates stellar M/L ratio as a function of 
color using functions of the form logio{M/L) = ax + bx x Color + C, where ax and bx are 
bandpass dependent constants and C is a constant that depends on the stellar initial mass 
function (IMF). For the V band Bell et al. (2003) find ax = -0.628 and 6a = 1-305; assuming 
a Kroupa (1993) IMF, they find C = —0.10. This yields an expression for the stellar mass 
in Mq for a given B — V color: 

M, = ^;;„^10-°-^28+i.305(B-F)-o.io^ ^3) 

where 

Vlum = 10-°-^(^-^-«2). (4) 

Here, vium is the luminosity parametrized in terms of absolute V magnitude. 

How reliable are stellar masses determined from this procedure? Clearly, the above 
relationship between M^ and B — V cannot apply to all galaxies, and must depend on the 
assumed stellar IMF, and range of ages, dust, and metallicity. However, it is encouraging to 
note that several studies (Bell et al. 2003; Drory et al. 2004; Salucci, Yegorova, & Drory 2008) 
find generally good agreement between masses based on broad-band colors and those from 
spectroscopic (e.g. Kauffmann et al. 2003) and dynamical (Drory et al. 2004) techniques. 
Typical errors are within a factor of two to three. Salucci et al. (2008) derive disk masses 
with both photometric and kinematic methods and find the two methods are equivalent on 
average. For a sample of 18, they find an rms scatter of 0.23 dex, while on an individual 
basis the deviation can be as high as 0.5 dex. 

We used this relation to compute stellar masses for 126 of 143 (88%) objects. The 
remainder did not have B — V colors available in the Hyperleda database or RC3. The mass 
distribution is summarized in Figure O Individual masses are listed in Table [TJ This sample 
of 126 galaxies is referenced henceforth as sample S2. 



3. Method and Analysis 

The structural properties of galaxy components, such as bulges, disks, and bars can be 
derived through the decomposition of the 2D light distribution, taking into account the PSF. 
There are several algorithms for 2D luminosity decomposition, including GIM2D (Simard et 



al. 2002), GALFIT (Peng et al. 2002), and BUDDA (de Souza et al. 2004). The latter two 
allow bulge-disk-bar decomposition, while the former only allows bulge-disk decomposition. 

Most previous work has addressed 2D bulge-disk decomposition only. Allen et al. (2006), 
for example, performed bulge-disk decomposition of 5-band images with GIM2D on 10,095 
galaxies from the Millennium Galaxy Catalog (Liske et al. 2003; Driver et al. 2005). However, 
recent work (e.g., Laurikainen et al. 2005; Graham & Balcells, in preparation) has shown 
that the B/T ratio can be artificially inflated in a barred galaxy unless the bar component 
is included in the 2D decomposition. The fact that most (> 60%) bright spiral galaxies are 
barred in the NIR (Eskridge et al. 2000; Laurikainen et al. 2004; MJ07; Menendez-Delmestre 
et al. 2007), further warrants the inclusion of the bar. Another advantage of bulge-disk-bar 
decomposition is that it allows us to constrain the properties of the bar itself, and to constrain 
scenarios of bar-driven evolution (see § [1]). 

Motivated by these considerations, several studies have tackled the problem of 2D bulge- 
disk-bar decomposition. Laurikainen et al. (2005, 2007) have developed a 2D multicompo- 
nent decomposition code designed to model bulges, disks, primary and secondary bars, and 
lenses; they apply Sersic functions to bulges and use either Sersic or Ferrers functions to 
describe bars and lenses. Reese et al. (2007) have written a non-parametric algorithm to 
model bars in ~ 70 /-band images. Gadotti & Kauffmann (2007) are performing 2D bulge- 
disk-bar and bulge-disk decomposition of 1000 barred and unbarred galaxies from SDSS with 
the BUDDA software. 

In this study, we perform 2D two-component bulge-disk decomposition and three- 
component bulge-disk-bar decomposition of the OSUBSGS sample with GALFIT. We note 
that Laurikainen et al. (2007) have also performed bulge-disk-bar decomposition on the OS- 
UBSGS sample. However, there are also important complementary differences between our 
study and theirs. The decomposition algorithm and tests on the robustness performed in our 
study are different (see §|3]and §|1]). Furthermore, unlike Laurikainen et al. (2007), we also 
compare the bulge-to-total ratio {B/T) with predictions from hierarchical models of galaxy 
evolution (§[5]), and also present the distribution of bar-to-total ratio (Bar/T). 



3.1. Image Preparation 

Running GALFIT on an image requires initial preparation. The desired fitting region 
and sky background must be known, and the PSF image, bad pixel mask (if needed), and 
pixel noise map must be generated. We addressed these issues as follows: (1) The GALFIT 
fitting region must be large enough to include the outer galaxy disk, as well as some sky 
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region. Since cutting out empty regions of sky can drastically reduce GALFIT run-time, a 
balance was sought between including the entire galaxy and some decent sky region, while 
excluding large extraneous blank sky areas. (2) It is possible for GALFIT to fit the sky 
background, but this is not recommended. When the sky is a free parameter, the wings of 
the bulge Sersic profile can become inappropriately extended, resulting in a Sersic index that 
is too high. Sky backgrounds were measured separately and designated as fixed parameters; 
(3) GALFIT requires a PSF image to correct for seeing effects. Statistics of many stars in 
each frame can be used to determine an average PSF. However, many of our images contain 
merely a few stars. Instead, a high S/N star from each frame was used as a PSF. (4) We used 
ordered lists of pixel coordinates to make bad pixel masks, which are useful for blocking out 
bright stars and other image artifacts. (5) We had GALFIT internally calculate pixel noise 
maps for an image from the noise associated with each pixel. Noise values are determined 
from image header information concerning gain, read noise, exposure time, and the number 
of combined exposures. 



3.2. Decomposition Steps 

Figure H] summarizes our method of decomposition, which we now detail. GALFIT 
requires initial guesses for each component it fits. It uses a Levenberg-Marquardt downhill- 
gradient algorithm to determine the minimum x^ based on the input guesses. GALFIT con- 
tinues iterating until the x^ changes by less than 5e-04 for five iterations (Peng et al. 2002). 
We recognize that a drawback to any least-squares method is that a local minimum, rather 
than a global minimum, in x^ space may be converged upon. We explore this possibil- 
ity with multiple tests described in § HI We adopted an iterative process, involving three 
separate invocations of GALFIT, to perform 1-component, 2-component, and 3-component 
decomposition: 

1. Stage 1 (single Sersic): In Stage 1, a single Sersic component is fitted to the galaxy. 
This serves the purpose of measuring the total luminosity, which is conserved in later 
Stages, and the centroid of the galaxy, which is invariant in later fits. 

2. Stage 2 (exponential plus Sersic): In Stage 2, the image is fit with the sum of an 
exponential disk and a Sersic component. During the Stage 2 fit, the disk h/a and PA 
are held constant at values, which we take from the published ellipse fits of MJ07, as 
well as ellipse fits of our own. This procedure reduces the number of free parameters 
in the fit by fixing the disk h/a and PA, which are easily measurable parameters. It 
also prevents GALFIT from confusing the disk and bar, and artificially stretching the 
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disk along the bar PA in an attempt to mimic the bar. As initial guesses for the 
Sersic component in Stage 2, the output of Stage 1 is used. The Sersic component in 
Stage 2 usually represents the bulge, in which case Stage 2 corresponds to a standard 
bulge-disk decomposition. 

However, in a few rare cases, where the galaxy has just a bar and a disk, the Sersic 
component in Stage 2 represents a bar. The latter is recognizable by a low Sersic index 
and large half-light radius. 

3. Stage 3 (exponential plus two Sersic components): In Stage 3, a three-component 
model consisting of an exponential disk, a Sersic bulge, and a Sersic bar is fit. As 
suggested by Peng et al. (2002), the bar can be well described by an elongated, low- 
index Sersic {n < 1) profile. As in Stage 2, the disk b/a and PA are held constant at 
values predetermined from ellipse fits. We provide initial guesses for the bar b/a and 
PA, based on ellipse fits of the images from MJ07 or analysis of the images in DS9. 
We provide GALFIT with input guesses for the bulge parameters, based on the output 
from Stage 2. In principle, it is also possible to generate reasonable guess parameters 
for the bulge and disk from a bulge-disk decomposition on a ID profile taken along 
a select PA. As described in § 14. 3[ we also experiment with initial guesses derived in 
this way, and find that the final convergence solution is the same. We also note that 
GALFIT fixes the bulge b/a and does not allow it to vary with radius, while real bulges 
may have a varying b/a. We tested the impact of fixed and varying bulge b/a on the 
derived B/T (§ 14.11) and find that there is no significant change in B/T. 

For objects with central point sources, the bulge Sersic index in the Stage 2 and Stage 
3 models can be inadvertently overestimated unless an extra nuclear component is added 
to the model. Balcells et al. (2003) note that for galaxies imaged both from the ground 
and HST, the combination of unresolved central sources and seeing effects mimic high- 
index bulge Sersic profiles in the ground images. Depending on sample and resolution, the 
frequency of central sources can range from 50% to 90% (Ravindranath et al. 2001; Boker 
et al. 2002; Balcells et al. 2003). Ravindranath et al. (2001) find a frequency of 50% in 
early type (E, SO, SO/a) galaxies, while Boker et al. (2002) measure a frequency of 75% for 
spirals with Hubble types Scd to Sm. Balcells et al. (2003) determine a frequency of 84% for 
SO-Sbc galaxies imaged with HST. Our dataset most closely resembles the latter sample, 
so we might expect that, as an upper limit, a similar fraction of our galaxies will need to be 
corrected with an extra compact component. 

We added point sources as third or fourth components to the initial models. For those 
cases where the model successfully converged with the extra component, the images were 
visually inspected to verify the presence of a central bright source. Sometimes, the model 
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converged to significantly different and unreasonable parameters for all components. Other 
times, the model would converge to a very dim point source without changing any of the other 
model parameters. Where new model parameters were not unreasonable and not identical 
to the case without the point source, the new model was adopted. This was the case for 
111 of 143 (78%) of our sample. The point sources contribute less than 1%, 3%, and 5% of 
the total luminosity 55%, 86%, and 95% of the time, respectively. As the point sources take 
up such a small fraction of the light distribution, their contribution is folded back into B /T 
in all cases where a point source was modeled. Inclusion of the point source reduced bulge 
index by ~ 0.8 for both barred and unbarred galaxies. Such a change is expected based 
on the above discussion. The decline in bulge index caused a minor decrease in B /T\ on 
the mean, this change was 1.04% for barred galaxies and 0.32% for unbarred galaxies. For 
barred galaxies, this light most primarily added to D /T rather than Bar/T. 

It is important to recognize the physical significance of the added nuclear components. 
We began by determining which objects show evidence for AGN activity. The sample was 
checked against the catalog of Ho et al. (1997), the Veron Catalog of Quasars & AGN, 12th 
Edition (Veron-Cetty & Veron 2006), and NED. Of the 111 objects fit with point sources, 43 
(39%) contain AGN. An additional 20 (18%) possess HII nuclei according to Ho et al. (1997) 
and visibly show bright compact nuclei. The remaining 48 (43%) probably contain neither 
AGN nor HII nuclei, but could house nuclear star clusters. For these objects, we visually 
inspected the images to ensure there was a bright compact source at the center. We are 
confident that the fitted point source components have physical counterparts in the data 
images. 

GALFIT also allows a diskiness/boxiness parameter to be added to any Sersic or expo- 
nential profile. We did not use this parameter for any bulge or disk profiles. Bars in general 
have boxy isophotes, and we could have included the diskiness/boxiness parameter in the bar 
profiles. However, it was found that adding boxiness to the bar did not change any model 
parameters, including fractional luminosities -B/T, D /T, and Bar/T, by more than a small 
percentage, even though the appearance of the residual images improved in some cases due 
to the change in bar shape. Consequently, we chose to neglect bar boxiness altogether. 



3.3. Choosing the Best Fit Between Stage 2 and Stage 3 

All objects in our sample were subjected to Stages 1, 2, and 3. Depending on whether 
a galaxy with a bulge is unbarred or barred, its best fit should be taken from the Stage 
2 bulge-disk decomposition or the Stage 3 bulge-disk-bar decomposition, respectively. For 
objects with prominent bars, it is obvious that the Stage 3 model provides the best fit. 
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However, it is more difficult to decide between Stage 2 versus Stage 3 fits in galaxies which 
host weak bars with no strong visual signature. In practice, we therefore applied the set of 
criteria below to each galaxy in order to select between the Stage 2 bulge-disk decomposition 
and Stage 3 bulge-disk-bar decomposition. Table [T] lists the model chosen for each galaxy. 
Table [2] summarizes the model parameters from the best fits. 

For completeness, we note that for the few rare galaxies (see § 13.21) that have just a bar 
and a disk, the choice of a final solution is between the Stage 2 bar-disk decomposition and 
Stage 3 bulge-disk-bar decomposition. The same guidelines below can be used to identify 
the best model. 

1. GALFIT calculates a x^ and x^ for each model. It was found that x^ almost universally 
declines between the Stage 2 and Stage 3 fits for a given object. This is because in the 
Stage 3 fit, five extra free parameters (bar luminosity, rg, Sersic index, 6/a, and PA) 
are added with the Sersic bar component, allowing GALFIT to almost always make 
a lower x^ model during Stage 3. However, this does not necessarily mean that the 
solution in Stage 3 is more correct physically. Thus, an increasing x^ w^^s interpreted 
as a sign that the Stage 3 fit should not be adopted, but a decreasing x^ was not 
considered as a sufficient condition to adopt Stage 3. 

2. In cases with prominent bars, a symmetric light distribution due to unsubtracted bar 
light was often found in the Stage 1 and Stage 2 bulge-disk residuals. This was strong 
evidence that the Stage 3 bulge-disk-bar fit be selected. NGC 4643 is shown in Fig- 
ure [5] because it has a particularly striking bar residual; the corresponding fit param- 
eters appear in Table [31 Note that in all figures and tables, we adopt the convention 
that PA values are positive/negative if they are measured from North counterclock- 
wise/clockwise. 

3. The Stage 2 and Stage 3 models were selected only so long as the model parameters 
were all well behaved. In unbarred galaxies, the Stage 3 model parameters might be 
unphysically large or small, in which case the Stage 2 fit was favored. Conversely, in 
galaxies with prominent bars, the bulge component of the Stage 2 bulge-disk fit tends 
to grow too extended in size. Addition of a bar in the Stage 3 bulge-disk-bar fit removes 
this artifact, giving a more physical solution. An extreme example of this situation is 
the barred galaxy NGC 4548, which has a prominent bar and a faint disk. The Stage 2 
fit, based on a Sersic bulge and exponential disk, is highly inadequate to describe the 
bulge, disk, and the bar because it leads to an extremely extended bulge. The Stage 
3 bulge-disk-bar fit, however, yields a believable fit with a prominent bar. The results 
of Stage 1, Stage 2, and Stage 3 are displayed in Figure Eland Table JH 
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4. Not all barred galaxies had unphysical Stage 2 models. Instead, the bulge could be 
stretched along the PA of the bar, giving the bulge a lower Sersic index and larger 
effective radius. A Stage 3 model that returned the bulge to a size and shape more 
representative of the input image was favored over the Stage 2 fit. Figure [7] and Table O 
demonstrate this behavior in NGC 4902. We distinguish this effect from cases like NGC 
4548 (Figure [6] and Table H]) where the Stage 2 fit is completely wrong. 

5. In cases where there was no bar, GALFIT can sometimes be enticed into fitting a bar 
to any existing spiral arms, rings, or the clumpy disks of late-type spirals. Stage 3 fits 
in these cases could be discarded by noting the resulting discrepancies in appearance 
between the galaxy images and the Stage 3 model images. Examples of false bars are 
shown in Figure [HI 

After fitting the whole sample and picking the best fit from either the Stage 2 bulge-disk 
decomposition or the Stage 3 bulge-disk-bar decomposition, we also performed the following 
extra tests. For our sample SI of 143 bright {Mb < —19.3) low-to-moderately inclined 
{i < 70°) spirals in the OSUBSGS survey, we determine the fraction (75 of 143 or ~ 52%) 
of spiral galaxies where a bulge-disk-bar decomposition was picked as the best fit for the H- 
band image. There are also eight galaxies with pure bar-disk fits. The H-hand bar fraction, 
which is defined as the fraction of disk galaxies that are barred, is therefore 58.0 ± 4.13% 
(83 of 143). We then compared our results (58.0% ± 4.13%) with the H-hand bar fraction 
(60%) determined from ellipse fits of the OSUBSGS sample by MJ07, with a slightly more 
conservative inclination cut {i < 60°). The two numbers are in excellent agreement. As 
a further check to our fits, we compare the bar and unbarred classification for individual 
galaxies from our fits with those from MJ07, which were based on ellipse fits. Of the 73 
galaxies that we classify as barred, and that are mutually fitted by MJ07, 54 (74%) are also 
classified as barred by MJ07. The remaining 19 (26%) galaxies are mainly weakly barred 
(with Bar/T below 0.08). Their RC3 optical types are weakly barred AB (10), barred B (7), 
and unbarred A (2). 

In most previous bulge-disk and bulge-disk-bar 2D decomposition, the issue of param- 
eter coupling and the systematic exploration of local versus global minima in x^ have been 
ignored. Quantifying how the parameters are coupled is important in measuring error bars 
for the model parameters. With 2D models containing several free parameters, this is not 
an easy task. Although we also do not address this problem in rigorous detail, we describe 
in § 14.41 simple test that explores parameter coupling in our models. 
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4. Extra Tests to Verify Correctness of Fits 

4.1. Varying h/a as a Function of Radius 

Models generated with GALFIT do not allow the h/a of the bulge, disk, or bar to vary 
with radius. Since real bulges may have a varying 6/a, it is legitimate to investigate the 
impact of fixing the bulge h/a on the estimated B /T. We therefore performed the following 
test on NGC 4548. To mimic a model bulge of varying h/a, we fitted the bulge light of 
NGC 4548 with ten concentric Sersic profiles of increasing re and varying h/a. The rg of 
the outermost profile comes from the original bulge model (see Table H]) where h/a was kept 
constant with radius. The separation in rg between adjacent profiles is 0.5 pixels (0.75"). 
The luminosity, Sersic index, h/a, and PA of each profile were free parameters. The disk 
and bar components were fixed to the values in Table HI as the emphasis was on the change 
in the bulge. 

Figure [9] compares the B /T obtained by fitting the bulge of NGC 4548 with a Sersic 
model of constant h/a as opposed to a Sersic model varying h/a. The bulge h/a (0.88), PA 
(-66.5), and B/T (13%) from the original Sersic fit of constant h/a (Table H]) are indicated 
with horizontal lines on the 3 panels. The top two panels show the run of h/a (0.85 to 1.0) 
and PA (—90° to +90°) of the ten concentric Sersic profiles. It can be seen that the Sersic 
indices of the ten bulge models were generally higher toward the center and declined at larger 
Te, indicating that the 'fitted bulge is more concentrated at the center. The bottom panel 
shows the cumulative B/T calculated by summing all models with r < r^. The last point 
representing the summed B/T from all ten components is 14.5%, in good agreement with 
the 13.0% value from the Sersic fit of constant h/a. Thus, using a Sersic model of constant 
h/a, does not have any significant adverse impact on our derived B/T in NGC 4548. 



4.2. Fitting Artificially Simulated Images 

An elementary test is to determine if GALFIT can recover the known parameters of 
artificially simulated noisy images. The images were simulated by taking parametric model 
images produced by GALFIT, and adding noise to the images with the PyFITS module 
for Python (Barrett & Bridgman 1999). Noise in each pixel was calculated by adding in 
quadrature the noise due to the source, sky, and read noise. The standard deviation of pixel 
noise in electrons was computed as 



<7 — yTsource + Tg^y + T^g^^, (5) 
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where Tgource is the number of electrons due to the source, Tsky is the number of electrons 
due to the sky, and Tread is the detector read noise. The contribution due to detector dark 
current was very small and therefore neglected. The offset added to each pixel was drawn 
from a normal distribution centered at zero with standard deviation a. 

Our test sample consisted of 40 models (20 bulge-disk and 20 bulge-disk-bar) with 
noise added as described. Thirty of the images included point sources as extra components. 
The range explored for each parameter in the model images is a fair representation of the 
parameter space covered by our full sample (e.g., B /T ranges from 0.02 — 0.70, the bulge 
index ranges from ~ — 5, and the full range of possible bulge and bar PA was also tested 
(—90° to -|-90°)). In terms of surface brightness, the models span five magnitudes in mean 
surface brightness inside the disk scalelength. Examples of the noise-added models are shown 
in Figure [TUl 

The noisy images were subjected to the 2D decomposition procedure outlined in Fig- 
ure HI GALFIT reproduced the model (bulge, disk, bar, and point source) parameters quite 
closely for the majority of the test cases. Figure [TT] compares the recovered versus original 
model parameters. Except for some extreme cases where the images were highly distorted 
by noise, all parameters were recovered to within a few percent. Figure [12] plots the ratio of 
model-to-recovered parameter against mean surface brightness inside the disk scalelength; 
there is no strong trend in error with dimming surface brightness. This suggests our de- 
compositions are effective across the parameter space spanned by our sample. The overall 
success of this test is evidence that GALFIT is able to converge to the absolute minimum in 
X^ space for our bulge-disk and bulge-disk-bar decompositions when the input is the sum of 
parametric functions. 



4.3. Using ID Decomposition To Generate Guesses for Bulge Parameters 

It is important to verify that GALFIT converges to the same solution even if the initial 
guesses for the bulge parameters in Stage 2 and 3 are different. Bulge-disk decomposition 
from ID profiles provides an alternative means of generating initial guesses. While ID 
bulge-disk decompositions of radial profiles along the bar major axis can be influenced by 
the bar, decomposition of cuts along the bar minor axis will not be influenced as heavily. 
The resulting bulge and disk parameters should be adequate guesses for Stage 3 of our 2D 
decomposition method. 

We tested the robustness of our Stage 3 fits by extracting initial guesses for the bulge 
and disk using ID decomposition along the bar minor axis. The nonlinear least-squares 
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algorithin designed to perform the ID decomposition simultaneously fits the sky-subtracted 
profiles with the sum of a Sersic bulge and an exponential disk, while ignoring the PSF. 
The results from the ID decomposition include a bulge magnitude, r^, Sersic index, disk 
magnitude, and disk scalelength. 

The robustness of several bulge-disk-bar fits were tested by using the results of the ID 
decomposition as input to Stage 3. The ID decompositions do not provide information about 
the axis ratio (b/a) or PA, so these parameters for the bulge were estimated by eye; for the 
disk, the b/a and PA were fixed to the values determined by ellipse fitting, as described in 
§ 13 .21 The initial bar parameters were unchanged from the earlier Stage 3 fits. In all cases, 
the new models were identical to the Stage 3 models. As an example. Table |6] compares 
Stage 3 input derived from ID decomposition and GALFIT for NGC 4548 and NGC 4643. 
In each case, both sets of input reproduced the same results. 



4.4. Parameter Coupling 

Assessing the coupling between model parameters is complicated when models have a 
large number of free parameters. A standard approach is to calculate confidence regions using 
multi-dimensional ellipsoids for a given Ax^ contour. As the errors in the GALFIT models 
are not normally distributed, but are instead dominated by the systematics of galaxy struc- 
ture, this approach does not yield meaningful results because of the ambiguity in assessing 
which Ax^ contour levels are statistically significant. 

We carry out a simple test for representative galaxies to determine not only the effects 
of parameter coupling, but also the the effect parameter coupling has on model parameter 
errors, paying particularly close attention to B/T, D/T, and Bar/T, as they are of primary 
interest. We perform this test on four representative galaxies (NGC 3885, NGC 4151, NGC 
4643, and NGC 7213). Two are barred (NGC 4151 and NGC 4643). Two have high {n > 2) 
bulge indices (NGC 4643 and NGC 7213), and the other two have low {n < 2) bulge indices. 

We fit two and three-component bulge-disk and bulge-disk-bar models using fixed bulge 
indices of n = 1 and n = 4. The initial inputs to these fits were the same as those used to 
generate the model in which bulge n is a free parameter. We then had GALFIT re-fit the 
models with bulge index as a free parameter using these output model parameters as input 
initial guesses. We compare the output of these two fits with the the best fit as selected in 
§ 13. 3[ in which bulge n is a free parameter. Ideally, the re-fits should converge to the same 
final parameters as the selected best fit. 

Table [7] displays the outcome of this test for the four representative systems. For barred 
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galaxies NGC 4643 and NGC 4151, the n = 1 and n = 4 re-fits converged to the same x^ 
as when bulge n is initially left a free parameter. B/T, D/T, and Bar/T are precisely equal 
for NGC 4643, while for NGC 4151 there is a small dispersion of 0.1-0.2%. For unbarred 
galaxies NGC 3885 and NGC 7213, the n = 1 re-fits again converged to the same x^ ^^^d 
model parameters as when bulge n is initially left a free parameter. This is not true for the 
n = 4 re-fits. When the n = 4 condition is enforced, the bulges in these cases become too 
extended and luminous while the x^ drop below those when bulge n is initially left a free 
parameter. During the n = 4 re-fits, B/T increases further at the expense of D/T and the 
X^ remain unchanged or decrease further. The B/D ratio from the n = 4 re-fit for NGC 
3885 is 1.3, roughly 3.5 times higher than when n is free. Given that NGC 3885 is an SO/a 
galaxy with a bulge embedded in a smooth extended disk, the latter B/D is arguably too 
large. For NGC 7213, the n = 4 re-fit yields starkly unphysical values. The B/D ratio is 
11.3, and r^/h, the ratio between bulge effective radius and disk scalelength is 9.0. The lower 
x'^ in these cases cannot be taken seriously as the bulges are too luminous and the resulting 
B/D ratios do not match the data images. 

As illustrated by the above discussion and Table [TJ this test shows that in some cases 
(e.g., NGC 4643 and NGC 4151) GALFIT converges to the similar model parameters and 
B/T, D/T, and Bar/T while starting from highly different initial input guesses (e.g. bulge 
n = 1, bulge n = 4, and bulge n based on the Stage 2 fits) in different regions of parameter 
space. For NGC 3885 and NGC 7213, however, during the n = 4 re-fits GALFIT converged 
to models that were unphysical and different compared with the reasonable models generated 
with input guesses corresponding to the bulge fixed at n = 1 or the bulge n based on output 
from Stage 1. In effect, when the initial input guesses were very different from the data 
images, the resulting models were found, in spite of the lower x^, to be unphysical through 
comparison with the input data images. We emphasize that for all sources analyzed in the 
paper, the data, converged model output, and residuals are always inspected before adopting 
the best final fit (see § 13.31) . 

Table [7] also provides hints as to how the model parameters are coupled to bulge index. 
As suggested in the above discussion, fixing the bulge index to n = 4 leads to a more 
extended and luminous bulge, causing bulge r^ and B/T to rise without fail for increasing 
bulge index. The disk is coupled with the bulge such that increasing bulge index, bulge rg, 
and B/T yields a reduction in D/T. At the same time, disk scalelength either increases 
(NGC 3885 and NGC 4643) or decreases (NGC 3885 and NGC 4151); in the latter two 
cases, the disk becomes very compact and the bulge quite extended. The behavior of the 
bar is coupled with both the bulge and disk. In the case of NGC 4643, as bulge index is 
raised to n = 4, bar r^ becomes slightly larger, but Bar/T falls by a factor of 2.6 as light is 
redistributed from the bar and disk to the bulge. Bar index also declines as the bar assumes 
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a flatter profile. For NGC 4151, bar rg again increases slightly, but this time Bar/T rises by 
a factor of 1.7 as light is transferred from the disk to the bulge and bar. 

Based on the above test, we stress that GALFIT was able to overcome this parameter 
coupling in the cases where the input guess parameters well-reflected the data images. 



5. Results and Discussion 

5.1. Impact of Bars in 2D Decomposition 

From the Stage 2 bulge-disk decomposition and Stage 3 bulge-disk-bar decompositions, 
which we performed on all objects (§ 13. 2p we saw firsthand the effects of adding a bar to 
the fit of a barred galaxy. We summarize below some of these effects in order to underscore 
the importance of including a bar component in the 2D luminosity decomposition of barred 
galaxies 

1. During the Stage 2 bulge-disk decomposition of a barred galaxy, the luminosity which 
comes from the galaxy's disk, bulge, and bar gets distributed only between two model 
components: the model bulge and disk. Since the disk h/a and PA are measured 
independently and held constant during the fits, the Stage 2 model tends to distort the 
bulge in order to fit the bar. Thus, the bulge in the Stage 2 bulge-disk decomposition of 
a barred galaxy can be artificially long or too bright and extended. When a model bar 
component is added in the Stage 3 bulge-disk-bar decomposition of a barred galaxy, 
it forces a reshuffling of the luminosity between the three components. Generally, the 
bulge declines in luminosity, whereas light can be either taken from, or added back, to 
the disk. 

2. We find that the inclusion of a bar component in the Stage 3 bulge-disk-bar decom- 
position of a barred galaxy reduces the bulge fractional luminosity B /T, compared 
with the Stage 2 bulge-disk decomposition. For our 75 barred galaxies, the reductions 
correspond to factors of less than two, 2 to 4, and above 4, in 36%, 25%, and 39% 
of barred galaxies, respectively. The larger changes in B/T occur in very strongly 
barred galaxies, where a prominent bar cause the Stage 2 bulge-disk decomposition to 
overestimate the bulge. For instance, B/T declines in both of NGC 4643 (Figure [5] and 
Table [3]) and NGC 4548 (Figure E] and Table H]). In the latter case, B/T is reduced by a 
factor of ~ 5 between Stage 2 and Stage 3. These examples underscore the importance 
of including bars in 2D luminosity decomposition of very strongly barred galaxies. 
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3. The scalelength of the disk is generally unchanged by including the bar. NGC 4548 
(Figure [6] and Table Hj) is a good example. Sometimes, however, the disk from the 
Stage 2 bulge-disk decomposition of a barred galaxy is erroneous due to a poor fit. 
The disk parameters from the Stage 3 bulge-disk-bar decomposition are quite different 
in such cases. NGC 4643 (Figure [5] and Table [3]) illustrates this behavior. 

We find that for our sample SI of 143 bright {Mb < —19.3) low-to-moderately inclined 
{i < 70°) spirals (Figure [I]) in the OSUBSGS survey, 75 of 143 or ~ 52% are better fit with 
a Stage 3 bulge-disk-bar decomposition than a Stage 2 bulge-disk decomposition. There are 
also eight galaxies with pure bar-disk fits. As stated in § 13. 3[ the resulting if -band bar 
fraction (58.0% ±4.13% or 83 of 143) is in excellent agreement with the H-hsmd bar fraction 
of 60% reported by MJ07 based on ellipse fits of the OSUBSGS sample, with a slightly more 
conservative inclination cut (i < 60°). 



5.2. Mass in Bulges, Disks, and Bars 

The fractional if -band luminosities in the bulge, disk, and bar {B/T, D/T, Bar/T) of 
each galaxy can be considered as a fractional mass if we assume that the same mass-to-light 
{M/L) ratio can be used to convert the if -band luminosities of both the numerator {B, D, 
or Bar) and the denominator (T) terms into a stellar mass. This is not an unreasonable 
assumption as the if -band M/L ratio is not very sensitive to differences in dust or age that 
might exist between the bulge, disk, and bar. The uncertainties in M/L can be estimated 
by looking at population synthesis models. Chariot, Worthey, & Bressan (1996) find that 
for idealized galaxies with a single generation of stars, the uncertainties in M/L ratio due 
to different input stellar models and spectra are roughly ±35% for a fixed metallicity and 
IMF. Furthermore, as the age of a stellar population varies from ~ 0.5 Gyr to 10 Gyr, the 
if -band M/L ratio rises by a factor of ~ 2 to 3 (Chariot 1996). Asymptotic giant branch 
(AGB) stars dominate the NIR light for ages between 0.1 and 1 Gyr, while red giant branch 
(RGB) and super-giant branch (SGB) stars dominate between 1 Gyr and 10 Gyr. 

In this paper, we convert the B/T light ratio determined from if -band images to a B/T 
mass ratio by assuming a constant mass-to-light {M/L) in the if -band for both the bulge and 
the rest of the galaxy. However, Moorthy & Holtzman (2006) present line strengths of bulges 
and inner disks for 38 spirals with Hubble types SO to Sc. They show 76% of spirals have 
negative metallicity gradients. B — K color gradients are shown to largely match metallicity 
gradients and are likewise negative outward, indicating bulge M/L is higher than in the disk 
and bar. If bulges are much older than the disks, then our prescription would underestimate 
the true B/T mass ratio. If we assume an extreme case where bulges are ~ 12 Gyr and the 
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disk light is dominated by a young 3 Gyr population, our assumption of a constant i7-band 
M/ L ratio would underestimate the true B/T by a factor of < 2 (Chariot 1996). In several 
sections of the paper (e.g., § 15.31 § 15. 8p . we illustrate how our main results would change 
if the true B/T was higher by up to a factor of two. On the other hand, central regions 
of galaxies may harbor intense episodes of star formation. If the bulge is younger than the 
disk and happens to have star formation and a significant young population of massive stars, 
then our prescription could overestimate the true B/T mass ratio. This would make our 
current results on the high fraction of low B/T bulges even stronger. 



Using the total galaxy stellar mass from § 12.21 the fractional masses can be converted 
into absolute masses. (We do not convert the if -band luminosity directly into a mass as 
the if -band images do not have photometric calibration). The results are shown in Table [Hi 
For our sample SI of 143 bright {Mb < —19.3) low-to-moderately inclined [i < 70°) spirals 
with a mass-weighted mean Hubble type of Sab-Sb, we find that ~ 70% of the stellar mass 
is in disks, ~ 10% is in stellar bars and ~ 20% is in bulges (with ~ 11% in n > 2 bulges 
and ~ 9% in n < 2 bulges). Thus while bulges with n < 2 are highly ubiquitous (see next 
section), they only account for a small fraction of the total stellar mass budget. 

Figure [T3l shows the stellar mass for bulges, disks, and bars along the Hubble sequence. 

It is useful to compare our results with those of Driver et al. (2006), who performed 
bulge-disk decomposition of i?-band images with GIM2D on 10,095 galaxies from the Mil- 
lennium Galaxy Catalog (Liske et al. 2003; Driver et al. 2005). They found 68.6% of the 
stellar mass to be in disks, and 32.6% in bulges (with 30.8% in high n bulges, and 1.8% in 
low n <2 bulges). Their study thus finds a higher stellar mass fraction in all bulges (32.6% 
vs our 18.9%), and in high n bulges (30.8% vs our 10.4%), and a lower fraction in low n <2 
bulges (1.8% vs our 8.4%), and disks-|-bars (68.6% vs our 71.6% + 9.6%). This difference can 
be attributed to the fact that the Driver et al. (2006) study did not perform bulge-disk-bar 
fits and thus, their B/T ratios may be skewed to higher values. 



5.3. Distribution of Bulge Index and B/T 

Figure [H] shows the individual and mean B/T and bulge Sersic index, plotted, as a 
function of Hubble type and galaxy stellar mass. Barred and unbarred galaxies are shown 
separately. Figure [TBI shows the relationship between bulge index and B/T. 

We first consider the B/T values in Figure [HI The mean B/T in barred galaxies is 
lower than in unbarred galaxies, but there is a large overlap in the individual values. The 
offset in the mean B/T of barred and unbarred galaxies reported here, agrees with the result 
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of Laurikainen et al. (2007; see § 15.41) on the same sample. We also note that B /T does not 
correlate with Bar/T (Fig. [T6l) : aside from the six galaxies with large Bar/T (> 0.3), most 
galaxies have low-to-moderate Bar/T and a wide range of B /T is seen at each Bar/T. This 
is reassuring and suggests that the bar fit is not arbitrarily biasing the B /T values. The 
distribution of Bar/T is further discussed in § 15. 5[ 

How does the B /T vary as a function of Hubble type and galaxy stellar mass? Bulges 
with very high B /T (> 0.4) exist primarily in galaxies with high mass (M^ > 6 x 10^° M©) 
and early types (SO/a to Sab). Bulges with very low B /T (< 0.1) lie primarily in lower mass 
galaxies with later morphologies (Sb to Sc). It is striking that ~ 69% of bright (Mb < —19.3^ 
low-to-moderately inclined (i < 70°^ spirals have B/T < O.2.- these bulges are pervasive and 
exist across the whole spectrum of SO/a to Sd. The results are summarized in Table O We 
shall return to this point in § 15. 8[ We note again that these B/T mass ratios were calculated 
assuming a constant M/L ratio in the H-hand for the bulge and disk components. As noted 
in § 15.21 if the bulge in these high mass spirals is much younger (older) than the disk and 
bar, then the B/T can be overestimated (underestimated) by up to a factor of two, and the 
limiting value of 0.2 for the B/T cited in the above fraction, would have to be modified in 
the extreme case to 0.1 and 0.4, respectively. 

Some of the low B/T < 0.2 values for six barred SO/a and Sa galaxies in Figure UM may 
at first look suspicious. Balcells et al. (2007) report the mean B/T for SO galaxies to be 
0.25, so much smaller B/T are potentially worrisome. OSUBSGS if-band images of these 
objects in Figure [T71 show a smooth extended disk around the bulge. It should be noted that 
Hubble types were originally assigned on a combination of criteria including disk smoothness 
and spiral arm topology in addition to the prominence of the bulge. It is likely that these 
galaxies were assigned early Hubble types due to their smooth extended disks, in spite of 
their low bulge-to-disk ratio. 

Similarly, some of the high B/T ~ 0.4 bulges in three of the Sc galaxies may at first 
seem odd. However, again, visual inspection of their image (Figure [T7|) reveals prominent 
spiral arms and dumpiness, which may explain why they were assigned late Hubble types. 

How does the bulge Sersic index n vary as a function of Hubble type, and galaxy stellar 
mass (Figure [H]), as well as B/T (Figure [15])? The results are summarized in Tabled Only 
a small fraction (~ 1%) of bright spirals have high n > 4 bulges; such bulges lie primarily in 
so/a to Sab, and have a large B/T > 0.2. A moderate fraction (~ 22%) have intermediate 
2 < n < 4 bulges; these exist in barred and unbarred SO/a to Sd, and their B/T spans a 
wide range (0.03 to 0.5) with a mean of 0.29. A strikingly large fraction (^~ 76%j of bright 
spirals have low n < 2 bulges; such bulges exist in barred and unbarred galaxies across a wide 
range of Hubble types, and their B/T varies from 0.01 to O.4, with most having B/T < 0.2. 
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5.4. Comparison With Other Work 

As an independent check of our decomposition method, we compare our results with 
independently published decompositions. 

We find our mean if -band B/ D (Figure [TS]) ratios are comparable to the i^-band B/ D 
derived with the ID bulge-disk decompositions of Graham (2001) and Trujillo et al. (2002). 
Like Graham (2001) and Trujillo et al. (2002), we find B /D is widely variable with Hubble 
type and that mean B / D steadily declines from Sa through Scd galaxies. Graham (2001) 
and Trujillo et al. (2002) find bulge indices are widely variable within a Hubble type, but 
they are in general > 1 for early types and < 1 for late types. We likewise find wide scatter in 
bulge index with n <1 bulges existing in both early and late types. Figure [18] is also in good 
agreement with the more recent results of Graham & Worley (2008), who have calculated, 
as a function of morphology, inclination and dust-corrected B/ D and bulge Sersic indices. 
They find B / D values are typically < 1/3. 

Another meaningful comparison can be made with Laurikainen et al. (2007) who, using 
their own 2D decomposition code, fit a hybrid sample containing some OSUBSGS galaxies 
combined with additional SO galaxies. One difference between their work and ours is that 
they typically model bars with a Ferrers function, but may sometimes use a Sersic profile, 
while we use only the latter. Also, they include additional components to model secondary 
bars or inner disks. They report a distinct offset in the mean B/T between barred and 
unbarred galaxies, which we confirm in Figure [TH Their mean B/T are similar to ours, and 
they conclude that pseudobulges exist throughout the Hubble sequence. The Sersic indices 
derived by Laurikainen et al. (2007) are likewise similar, on the mean, to ours for both barred 
and unbarred systems. They likewise find n < 2 bulges across early and late Hubble type 
galaxies. 

Balcells et al. (2003) emphasized that bulges typically have indices ~ 3 or lower. Our 
results in Figure O are consistent. We find a low frequency (~ 1%) of high n > 4 bulges, 
with most bulges having n < 3. 



5.5. Bar Strength 

Stellar bars exert gravitational torques on the gaseous component and are particularly 
efficient in driving gas from the outer disk into the inner kiloparsec (see § [1]). Thus, it 
would be natural to have a measure of bar strength that is sensitive to the strength of the 
gravitational torque and hence measures both the shape and mass of the bar. 
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Many measures of bar strength have been formulated. The Qb method of (Block et 
al. 2002; Buta et al. 2003; Buta et al. 2005) measures directly the gravitational torque at a 
single point along the bar. This method requires a value of scaleheight for the disk and a 
model of the potential to be made from the image. In the bar /inter-bar contrast method of 
Elmegreen & Elmegreen (1985) and Elmegreen et al. (1996), bar strength is parameterized as 
the ratio between peak surface brightness in the bar region and the minimum surface bright- 
ness in the inter-bar region. Elmegreen & Elmegreen (1985) and Elmegreen et al. (1996) also 
characterize bar strength with the maximum amplitude of the m = 2 mode from Fourier 
decomposition. When ellipse fitting is apphed, the maximum ellipticity of the bar, Cbar, can 
be used to characterize bar strength (e.g. MJ07). This constitutes a partial measure of bar 
strength only, however, as it offers no information about mass of the bar. 

Bulge-disk-bar decomposition in the if-band provides another partial measure of bar 
strength through the if -band Bar/T light ratio, which is a measure of the Bar/T mass ratio 
under the assumption that the H-hand M/L ratio is the same for the bar and the rest of 
the galaxy, as discussed in § 15.21 Figures [19] and [20] explore the derived bar properties. 

The upper left panel of Figure [H] plots the individual and mean Bar/T against Hubble 
type. There is a wide range (~ 0.03 to ~ 0.47) in the individual Bar/T at a given Hubble 
type. The mean Bar/T remains fairly constant with Hubble type from Sa to Sb, but shows 
a possible weak decline by about 0.1 from Sb to Sc. Their number statistics are too small to 
make any robust statement for later Hubble types. We also note that six systems have high 
Bar/T above 0.3; these are displayed in Figure [211 

Bar Sersic indices are mostly below unity. Neither the individual, nor the mean bar 
Sersic index, show any trend with Hubble type or with stellar mass, for Sa to Sc galaxies 
(Fig. [THj). Thus, the steepness of the bar profile does not seem to depend on the Hubble type. 
Is the bar mass ratio and its mass profile related? There is a wide range in the individual 
Bar/T at a given bar Sersic index (Fig. [20]) . The mean Bar/T rises with bar index out to 
a bar index of ~ 0.6, and then flattens out. This suggests that on the mean, bars of lower 
Bar/T have flatter profiles. 

Is there a relation between the bar strength and the bulge present in a galaxy? There 
is a wide range in the individual Bar/T at a given B/T, and at a given bulge Sersic index 
(Fig. [20]) . The mean Bar/T shows a weak decline for bulge Sersic indices above two. Similarly 
the mean Bar/T shows a weak rise from 0.1 to 0.25 as B/T rises out to 0.15, after which 
the trend flattens or reverses. 

Both Bar/T and maximum bar ellipticity Cbar are partial measures of bar strength. 
Figure [19] shows mean Bar/T may scale weakly with Hubble type. The bars with highest 
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Cbar (i-e, thin bars) are often termed strong bars, and Cbar has been shown to correlate with 
Qb- Total bar strength should scale with both bar mass and bar ellipticity. Does bar strength 
have a dependence on Hubble type? The upper left panel of Figure [20] plots the product of 
Bar/T and Cbar, as determined by MJ07 for galaxies mutually classified as barred, against 
Hubble type. There is a wide range in Bar/T x Cbar in each bin, and mean bar strength 
shows no definite trend with Hubble type. We note that bars with high Bar/T and high Cbar 
should exert the largest gravitational torque and be most effective at driving gas inflows. A 
nice example is the oval or lens galaxy NGC 1317 (Figure [2T]) : the bar has a low ellipticity, 
but its B/T is large as it is extended and massive. Such bars or lenses may exert significant 
gravitational torques although they are not very elongated. 



5.6. Bar Fraction as a Function of B/T and Bulge Index 

As outlined in § 15. 3[ we found that as many as ~ 76% of bright spirals have bulges with 
n <2; such bulges exist in barred and unbarred galaxies across a wide range of Hubble types, 
and their spread in B/T is from 0.01 to 0.4, with most having B/T < 0.2. The variation 
of the bar fraction as a function of B/T and bulge n can provide important constraints on 
bulge formation scenarios (§ 15. 8p . Table [10] shows our results. The bar fraction declines 
with bulge index; ~ 65% of the spirals with low n < 2 bulges host bars while intermediate 
2 < n < 4 bulges have a lower bar fraction (~ 38%). The high n > 4 bulges in the sample 
are unbarred, so the bar fraction is 0%. Systems with low B/T are more likely to be barred. 
For B/T < 0.2, the bar fraction is high (~ 68%). Systems with 0.2 < B/T < 0.4 and 
B/T > 0.4 have lower bar fraction (~ 42% and ~ 17%). 

Overall, Table [TU] shows bulges with low n < 2 and low B/T < 0.2 preferentially exist 
in barred galaxies. This is consistent with earlier work (Odewahn 1996; Barazza et al. 2008; 
Marinova et al. 2008; Aguerri et al. 2008, in prep.) where an enhanced optical bar fraction is 
seen is galaxies with late Hubble types or low B/D. It may be tempting to infer this result 
to mean the formation pathway of two-thirds of low-B/T bulges is related to bars in that 
spontaneous or/and tidally induced bars play a role in bulge formation (with the remaining 
one-third of such bulges may have been formed either by mechanisms like retrograde minor 
mergers or short-lived bars). We caution that this type of cause-effect relationship is not 
the only scenario consistent with this result. It may also be possible that bar instabilities 
are favored in galaxies with low B/T and no inner Lindblad resonances (ILR). Under these 
conditions, the swing amplifier model with a feedback loop (Julian & Toomre 1966; Toomre 
1981; Binney & Tremaine 1987) may be responsible for bar formation and partly account 
for the high bar fraction in galaxies of low B/T. 
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5.7. Formation of Bulges 

Our observational results provide some interesting challenges for models of galaxy evolu- 
tion that try to address the origin of present-day bulges. Any successful model must be able to 
account for the observed distribution of bulge B/T and n in high mass (M^ > 1.0 x IO^^Mq) 
spirals, as shown in Table and Table [101 In particular, the following results must be 
reproduced: 

1. In terms of the overall distribution of bulge n, as much as (~ 74%) of high mass spirals 
have bulge n < 2: such bulges exist in barred and unbarred galaxies and their B/T 
ranges from 0.01 to 0.4, with most having B/T < 0.2 (Table [9]). A moderate fraction 
(~ 24%) of high mass spirals have 2 < n < 4, and just (~ 2%) have n > 4. 

2. Theoretical models often make more robust predictions on the bulge-to-total mass 
ratio B/T than on the bulge index n, so we consider the empirical B/T distribution 
in detail. We note that as much as ~ 66% of high mass spirals have bulges with 
B/T < 0.2 (Table [9]). In terms of bar fraction, ~ 68% are barred (Table [TOj) . 

3. The fraction of bars rises among spirals with low bulge index n. About 63% of spirals 
with low n <2 bulges host bars, while the bar fraction in spirals with 2 < n < 4 bulges 
(44%) is two-thirds as large (Table [TO!). 

In a hierarchical Universe, there are several physical processes that contribute to the 
assembly of bulges: major mergers, minor mergers, and secular evolution. We briefly describe 
these, expanding on our introduction in § [H 

Major mergers, defined as those with mass ratio M1/M2 > 1/4, typically destroy the 
extended outer stellar disks during violent relaxation, leaving behind a classical bulge. Such 
bulges are associated with modest-to-high bulge Sersic indices, in the range 2 < n < 6 
(Hopkins et al. 2008; Springel et al. 2005; Robertson et al. 2006; § l5.8p in simulations. This 
trend is also consistent with the fact that among ellipticals, high luminosity ones tend to 
have a Sersic index n > 4, while low luminosity ones tend to have 2 < n < 3.5 (Caon et 
al. 1993; Kormendy et al. 2008). The final Sersic index depends on the amount of residual 
gas the settles into a somewhat disky component. Simulations by Hopkins et al. (2008) 
find that the Sersic indices of remnants from 1:1 gas-rich major mergers lie in the range of 
2 < n < 4, with most above 2.5 (see Fig. [22]). This body of evidence strongly suggests that 
many bulges with n > 2 might have a major merger origin. 

Minor mergers, typically defined as those with mass ratio 1/10 < M1/M2 < 1/4, do 
not destroy the stellar disk of the primary system, but can contribute to building bulges via 
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three pathways. Firstly, a fraction Fsat of the sateUite's stellar mass can end up in the central 
region of the primary galaxy. The value of Fgat depends on how centrally concentrated the in- 
falling satellite is. Typically, the more diffuse outer stellar body is tidally stripped, while the 
central core sinks by dynamical friction to the central region (e.g., Quinn et al. 1993; Walker 
et al. 1996). Secondly, a non-axisymmetric feature (e.g., a stellar bar or bar-like feature) 
can be induced in the main disk, and gravitational torques exerted by the feature can drive 
gas into the inner kpc (e.g., Hernquist & Mihos 1995; Jogee 2006 and references therein), 
where subsequent SF forms a compact high v/a stellar component, or disky pseudobulge. 
Most of the gas inflow happens during the merger phase and large gas inflow rates (e.g., 
^ 1 Mq per year) may be generated. Thirdly, gas inflow can also be caused by direct tidal 
torques from the companion (e.g., Hernquist & Mihos 1995; Eliche- Moral et al. 2006). It is 
to be noted that in the simulations by Hernquist & Mihos (1995), the gas inflow generated by 
non-axisymmetric features (e.g., bar-like features) in the inner part of the disk is much larger 
than that caused by direct tidal torques from the satellite. In the recent work of Eliche- Moral 
et al. (2006), N-body simulations of minor mergers followed by flts of ID Sersic+exponential 
models to the remnants, suggest that the bulge Sersic index and B / D ratio can grow as a 
result of the central re-concentration of stellar disk material in the primary system by tidal 
forces. Minor mergers are frequent under ACDM, and the likelihood of multiple successive 
minor mergers occurring during the formation of a galaxy is high. Bournaud, Jog, & Combes 
(2007) study the effects of repeated minor mergers on galaxy structure. They show that a 
disk galaxy undergoing successive minor mergers will eventually transform into an elliptical 
galaxy with an r^/*^-law proflle and high V/a. However, galaxy growth is not completely 
merger-driven and the efficiency of minor mergers at creating ellipticals must be regulated 
by other mechanisms (e.g., cold gas accretion). 

In addition, the process of secular evolution can build a disky bulge (pseudobulge) 
between merger events. Here a stellar bar or globally oval structure in a non-interacting 
galaxy drives gas inflow into the inner kpc, where subsequent star formation forms a compact 
high v/a stellar component (e.g., Kormendy 1993; Jogee 1999; KK04; Jogee, Scoville, & 
Kenney 2005; Athanassoula 2005; Kormendy & Fisher 2005; Kormendy 2008). This process 
is different from that of minor mergers in the sense that it happens in the quiescent phase 
of the galaxy, between minor or major merger events. The prevalence of pseudobulges in 
galaxies of different Hubble types is discussed in KK04, and select examples of SO galaxies 
with pseudobulges are also shown in Kormendy & Cornell (2004) and KK04. 

The present-day bulge mass can be written as the sum of mass contributed from each 
process: 

^^^bulge -^^J-bulge ^ 1,/inaj ~r /mini ~r /min2 ~r /min3 ~r /sec J) \S^ ) 
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• /maj is the percentage of the bulge stellar mass, which is built by major mergers, 

• /mini is the percentage of the bulge stellar mass, which is built during minor mergers 
from stars of the satellite. This depends on the fraction Fgat of the satellite's stellar 
mass, which ends up in the central region of the primary galaxy during each minor 
merger. 

• /mm2 is the percentage of the bulge stellar mass, which is built during minor mergers 
from gas inflow caused by a tidally induced bar. 

• /mins is the percentage of the stellar mass, which is built during minor mergers from 
gas inflow caused by tidal torques from the companion. 

• /sec is the percentage of the stellar mass, which is built secularly between merger events 
from gas inflow caused by bars or ovals 

In § 15. 8[ we compare our derived distribution of bulge n and B /T with hierarchical 
models that model major and minor mergers, but not secular evolution. The main goal 
of the model is to see whether bulges built via major mergers can account for the large 
fraction of high mass spirals with bulges of low B /T or/and low n. A secondary goal is 
to see if a first order simplified prescription for minor mergers can broadly account for the 
observations. We stress here that bulge-building during minor mergers is modeled in a very 
simple way: all the stars in the satellite are assumed to contribute to the bulge of the larger 
galaxy (i.e., Fgat = 100%), and bulge-building via gas inflow driven through tidal torques 
and via gravitational torques from induced bars are ignored (i.e., /min2 = 0, and /mins = 0). 
Furthermore, the models entirely ignore secular evolution between mergers. In a future 
paper, these extra terms will be addressed and a comprehensive picture built of the relative 
importance of minor mergers and secular processes in making present-day bulges. 



5.8. Comparison of B /T With Hierarchical Models of Galaxy Evolution 

We compare our data with the predictions from cosmological semi-analytical models 
based on Khochfar & Burkert (2005) and Khochfar & Silk (2006). We briefly describe the 
models first. The merger trees of dark matter (DM) halos are derived by using the extended 
Press-Schechter formalism (Press & Schechter 1974), as in Somerville & Kolatt (1999). When 
two DM halos merge, the merger time scale of the galaxies is calculated by considering the 
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timescale it would take the satellite galaxies to reach the central galaxy at the center of the 
halo via dynamical friction (e.g., Kauffmann et al. 1999; Springel et al. 2001). The baryonic 
physics, which includes radiative cooling, star formation, and feedback from supernovae, is 
treated via semi-analytic prescriptions (see Khochfar & Silk (2006) and references therein). 
Baryonic mass inside the dark matter halos is divided between hot gas, cold gas, and stars. 
The hot gas is initially shock-heated to the halo virial temperature. As the gas radiatively 
cools, it settles down into a rotationally supported disk at the halo center. Cold disk gas 
is allowed to fragment and subsequently form stars according to the Schmidt-Kennicutt 
law (Kennicutt 1998). Star formation is regulated by feedback from supernovae using the 
prescription in Kauffmann et al. (1999). 

Major mergers are typically considered as those with stellar mass ratio M^jMi > 1/4. 
In the simulations, one assumes that during a major merger any existing stellar disk is 
destroyed, gas is converted to stars with some star formation efficiency (SFE), and all stars 
present undergo violent relaxation to form a bulge. Therefore, the bulge-to-total stellar mass 
ratio (B/T) of a bulge immediately after a major merger is always one. Note that the SFE 
during a major merger is not assumed to be 100% as there is mounting evidence from SPH 
simulations (Springel & Hernquist 2005; Cox et al. 2008) that not all cold gas is converted 
to stars. Instead, the burst efficiency defined by Cox et al. (2008) is applied to control the 
fraction of stars formed due to the interaction. This efficiency is dependent on the relative 
masses of merging galaxies and is expressed as 

/ M Satellite \ /„n 

\^ivi Primary / 

where ei:i is the burst efficiency for a 1:1 merger and 7 fixes the dependence on mass ratio; 
Cox et al. (2008) find ei:i = 0.55 and 7 = 0.69. The remaining fraction (1-e) of gas is added 
to the gaseous disk and can start making stars. 

As stated above, immediately after a major merger, the remnant is a bulge with a B/T 
equal to one. As time proceeds, B/T falls because a stellar disk grows around the bulge as 
hot gas in the halo cools, settles into a disk, and forms stars. The formation of stars by any 
residual cold gas left at the end of the major merger also helps to grow the disk. Thus B/T 
falls until the next major merger happens, at which point B/T is reset to one in the models. 

The bulge may also grow in stellar mass due to minor mergers. Minor mergers are 
defined as mergers with mass ratio 1/10 < M1/M2 < 1/4, and the stellar disk of the large 
companion is not destroyed during such mergers. The models assume that during minor 
mergers, all the stars in the satellite are added to the bulge of the host, while the gas settles 
in the disk. When DM halos grow by accretion or minor mergers, the hot gas that comes 
in with a satellite is immediately stripped and added to the hot gas component of the host. 
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The cold gas in the disk of the sateUite is only added to the cold gas of the host if they 
merge. Until they merge the satellite is using up its own cold gas to make stars. 

Fig. [23] shows the relationship between the present-day B/T of a a high mass (M^, > 
1.0 X IO^^Mq) spiral and the redshift 2;iast of its last major merger. As expected, systems 
where the last major merger occurred at earlier times, have had more time to grow a disk 
and have a lower B /T . The dispersion in the present-day B/T at a given ziast is due to the 
different times spent by a galaxy in terms of being a satellite versus a central galaxy in a DM 
halo, since the cooling of gas and the growth of a disk is stopped when a galaxy becomes a 
satellite. Thus, galaxies that became a satellite galaxy shortly after their last major merger 
stayed at high B/T. Conversely, galaxies that continued to be a central galaxy for a long 
time after their last major merger will have low B /T . 

The present-day B/T of a high mass (M^ > 1.0 x IO^^Mq) spiral depends on its major 
merger history. In particular, we note from Fig. [23] that a high mass (M^ > 1.0 x lO^^M©) 
galaxy, which has undergone a past major merger since z <2 will end up hosting a present- 
day B/T > 0.2. In effect, a high mass spiral can have a present-day B/T < 0.2 only if its 
last major merger occurred at z > 2 (lookback times > 10 Gyr). 

The predicted distribution of present-day B/T depends on the galaxy merger history in 
the models and it is relevant to ask how well the latter is constrained observationally. Over 
the redshift range z ~ 0.24 to 0.80 (lookback times of 3 to 7 Gyr), recent studies by Jogee et 
al. (2008, 2009) find that among high mass (M^ > 2.5 x lO^^M©) galaxies, ~ 10% of galaxies 
are undergoing mergers of mass ratio > 1/10, and ~ 3% are undergoing major mergers of 
mass ratio > 1/4. These findings agree within a factor of less than ~ 2 with the merger rates 
from the models of Khochfar & Burkert (2001) over z ~ 0.24 to 0.80. At higher redshifts, the 
empirical merger rate/fraction is uncertain due to relatively modest volumes and bandpass 
shifting effects, but there is a general trend towards higher merger fractions at higher redshifts 
(e.g., Conselice et al. 2003). The models used here (Khochfar & Burkert 2001) agree with 
this trend and predict that ~ 13.5% and ~ 20% of high mass (M^ > 1.0 x IO^^Mq) spirals 
have undergone major mergers since z <2 and z < 4, respectively (see Table [TTj) . 



The contribution of galaxies with different merger histories to the present-day B/T 
distribution are shown in Table [11] The top and middle parts of the table describe systems 
with and without major mergers since z <2 and z < 4, respectively. In the model, ~ 13.5% 
of present-day high mass (M^ > 1.0 x IO^'^Mq) spirals, experienced a major merger since 
z < 2, causing most of them (~ 11.2%) to have a present-day high B/T > 0.4 and a negligible 
fraction (~ 0.1%) to have a low present-day B/T < 0.2. In contrast, the remaining ~ 86.5% 
spirals experienced no major merger at 2; < 2, and most (67.2%) of them have a present-day 
low B/T < 0.2. If the comparisons are extended to systems without a major merger since 
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z < 4, the numbers are very similar (see middle part of Table [TT]) . 

Table [TT] shows that there is good agreement between the model and data for the fraction 
of high mass spirals with present-day low B/T < 0.2 (~ 67% in the model versus ~ 66% in 
the data). The model contribution to low B/T < 0.2 comes almost entirely from galaxies, 
which have not had a major merger since z <2 (see column 4 in Table [TTj) . In fact, most of 
these galaxies have not even had a major merger since 2; < 4, as illustrated by the bottom 
part of Table [TTJ In the model, the fraction f~ 1.6%; column 3 of Table [77]) of high mass 
spirals, which have undergone a major merger since z < 4 and host a bulge with a present- 
day B/T < 0.2, is a factor of over thirty smaller than the observed fraction (^ 66%j of 
high mass spirals with B/T < 0.2. Thus, bulges built via major mergers since z < A 
seriously fail to account for the bulges present in ~ 66% of high mass spirals. These results 
are also illustrated in Fig. [231 which shows the comparison between data and models for the 
cumulative fraction of high mass spirals as a function of present-day B/T. 

It is also interesting to note from Table [11] that although the models reproduce well 
the frequency of bulges with present-day low B/T < 0.2, they tend to over-produce the 
frequency of present-day high B/T > 0.4 systems by nearly a factor of two (~ 14% in the 
model versus ~ 8% in the data; see columns 6 and 2 in middle part of Table 11). Most of 
this overproduction stems from major mergers at ;z < 4 (~ 13%; column 4 in middle part of 
Table [TT]) . This suggests that major mergers, as currently modeled here, are building bulges 
too efficiently. 

One possible solution to this problem might relate to the suggestion by Hopkins et 
al. (2009) that the efficiency of bulge-building during the major merger of two spirals depends 
not only on the mass ratio M1/M2 of the merger, but also depends on the cold gas mass 
fraction /gas in the disk. In their semi-analytic models, the entire stellar mass of the satellite 
violently relaxes, but the fraction of stellar mass in the primary disk that violently relaxes 
and adds to the bulge is M1/M2. This differs with our models where the entire primary 
stellar disk is always destroyed in a major merger. Furthermore, the fraction F of the total 
gas mass, which loses angular momentum, falls to the nucleus, and is transformed into stellar 
mass in a nuclear starburst, is ~ (1 — /gas) x {M1/M2). In particular, F is lower for more gas- 
rich disks, causing a suppression of the burst efficiency in gas-rich systems, and a reduction 
in the stellar mass that ends up in the bulge built during the major merger. The predictions 
for the distribution of B/T from the Hopkins et al. (2009) models are shown in the bottom 
of Table [TT] with major mergers considered as those with baryonic galaxy-galaxy mass ratio 
> /1/3. Due to the reduced stellar mass that ends up in the bulge after a major merger, the 
Hopkins et al. (2009) models tend to yield lower B/T after such a merger. Thus, the models 
produce overall fewer high B/T > 0.4 systems, and more intermediate 0.2 < B/T < 0.4 and 
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low B/T < 0.2 systems (see column 6 in Table [TTj) . Nonetheless, the predictions are not 
very different from the models by Khochfar & Burkert (2005) and Khochfar & Silk (2006), 
which we use in this paper. In particular, the conclusion that the large frequency of high 
mass (M^ > 1.0 x IO^^Mq) spirals with low B/T < 0.2 bulges can only be accounted for by 
spirals without a major merger since z <2 also holds with the Hopkins et al. (2009) models. 

Thus, we conclude that the observed large frequency (^ 66%) of high mass (M^ > 
1.0 X IO^^Mq) spirals with low present-day B/T < 0.2 can be accounted for in our and other 
hierarchical models by high mass spirals, which have not undergone a major merger since 
z < 2, and most of which have not even experienced a major merger over the last 12 Gyr 
since z < 4. In the models, most of these present-day low B/T < 0.2 bulges are built by 
minor mergers since z < 4. As noted earlier, our models explore bulge-building via minor and 
major mergers, but do not explicitly incorporate secular processes (see § 15. 7p . In practice, 
secular processes may contribute to the building of present-day low B/T < 0.2 bulges, and 
are particularly relevant at z < 2, where major mergers cannot build such bulges. 

For completeness, we further explore how sensitive are our results to other assumptions 
made in the data and models of Khochfar & Burkert (2005) and Khochfar & Silk (2006): 

• How sensitive are the results to the mass ratio used to separate major and minor 
mergers? Fig. [25] is similar to Fig. [21] except that the model now defines major mergers 
as those with mass ratio M1/M2 > 1/6. In this case, about 30% of the model spirals 
undergo major mergers since z < 4 rather than ~ 20%. The overall model F (black 
dashed line) now under-predicts the data F by about 10% for B/T > 0.2. However, 
the main conclusion that bulges built by major mergers since z < 4 cannot account 
for most of the low B/T < 0.2 bulges, present in a large percentage (~ 66%) of spirals 
still holds. 

• How sensitive are the results to the B/T cut used to define spirals? Fig. [26] is similar 
to Fig. [21] except that here spirals are considered to be systems with a B/T < 0.55 
rather than 0.75 in the models, and a corresponding cut is applied to the data points. 
The results are similar to Fig. [21] 

• How sensitive are the results to our assumed constant iJ-band mass-to-light (M/L) for 
the bulge, disk, and bar? Fig. [27] is similar to Fig. [21] except that the B/T of all the 
observed galaxies has been multiplied by a factor of two, in order to test what would 
happen in the case where the M/L ratio of the bulge in i7-band is twice as high as that 
of the disk and bar. This could happen in an extreme example where the dominant 
bulge stellar population was much older (e.g. 12 Gyr) than the age of the dominant 
disk stellar population (e.g., 3 Gyr). In such a case, the fraction of high mass spirals 
with B/T < 0.2 would change from ~ 66% in Fig. [21] to ~ 55%, and deviate from the 
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model overall model F (black dashed line) by ~ 20%. However, the main conclusion 
that bulges built by major mergers since z < 4 cannot account for most of the low 
B/T < 0.2 bulges, present in a large percentage (~ 55%) of spirals still holds. 

Finally, it is important to note that so far we have compared the data and model only 
in terms of bulge -B/T, but not in terms of bulge index n or in terms of bar fraction. In 
effect, we have only shown that the models reproduce a subset of the results outlined in 
points (l)-(3) of § 15.71 Since the semi-analytic models do not produce a distribution of bulge 
index ra, we resort to presenting only an indirect comparison in Table [121 We assume that 
bulges, which form in major mergers have a bulge n > 2.5. This assumption is based the 
evidence presented in § 15.71 Thus, in Test 1 of Table fT2l we compare the fraction (~ 66%) of 
galaxies in the semi-analytic models having B/T < 0.2 and no major merger since z < 4, to 
the observed fraction (~ 65%) of galaxies with B/T < 0.2 and bulge n < 2.5. There is close 
agreement between the two values. In Test 2 of Table [T2| the fraction (~ 12.7%) of model 
galaxies with B/T > 0.4 and no major merger since 2; < 4 is a factor of ~ 3 higher than the 
fraction (~ 3.5%) of high mass spirals with B/T > 0.4 and bulge n > 2.5. Similarly, the 
fraction (~ 1.6%) of model galaxies with B/T < 0.2 and a past major merger since z < 4 
is also a factor of ~ 2 higher than the observed fraction (~ 0.9%) of high mass spirals with 
B/T < 0.2 and bulge n > 2.5 (Test 3, Table [T2|) . Thus, in terms of bulge B/T and n, there 
is good agreement between data and model for Test 1 (involving model galaxies with no 
major mergers since z < 4). This suggests that the vast majority of bulges with B/T < 0.2 
and n < 2.5 likely formed in galaxies having had no major merger since z < 4. However, the 
agreement is less good for Tests 2 and 3 (involving model galaxies with major mergers since 
z < 4) and this suggests that the models may be building bulges a little bit too efficiently 
during a major merger, in agreement with the conclusion reached earlier. 

What about the role of bars in the formation of these bulges of low B/T and low n? 
A detailed direct comparison with the semi-analytic models is not possible as the role of 
bars is not yet modeled, but related comparisons are possible. First, it is important to 
note that bar-driven gas inflow into the inner kpc and the subsequent building of disky 
stellar components or 'pseudobulges' (see § [T]) can happen in both isolated galaxies and 
in minor mergers (§ 15. 7p . since bars can be spontaneously induced in an isolated disk or 
tidally induced during an interaction or minor merger. The triggering of a bar is favored 
in a prograde interaction or minor merger. Thus, bulge-building via induced bars is more 
likely to happen in prograde rather than retrograde minor mergers. Statistically about half 
of minor mergers might be prograde or prograde-like, and half retrograde. Thus, one would 
expect bars to be induced in only half of the minor mergers. If this assumption is correct 
and if most of the mass in bulges with present-day B/T < 0.2 is formed in minor mergers, 
then one would expect only about half of these bulges to host bars. This is close to what is 
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observed, as shown by Table [TOl We see that ~ 63% of high mass spirals with low n < 2 
bulges and ~ 68% of spirals with low B /T < 0.2 bulges host bars. This suggests that in 
high mass spirals, spontaneous and/or tidally induced bars may play a part in forming up to 
two-thirds of B /T < 0.2 or n < 2 bulges. The remaining one-third of such bulges may have 
been formed either by mechanisms that do not involve bars (e.g., retrograde minor mergers) 
or by bars that are not long-lived. 



6. Summary 

The properties of galaxy components (bulges, disks, and bars) in the local Universe 
provide key constraints for models of galaxy evolution. Most previous 2D decompositions 
have focused on two-component bulge-disk decomposition, and ignored the contribution 
of the bar even in strongly barred galaxies. However, as shown by this work and other 
recent studies (e.g., Laurikainen et al. 2005; Laurikainen et al. 2007; Reese et al. 2007), it 
is important to include the bar component in the 2D decomposition, in order to correctly 
estimate the bulge-to-total ratio (B/T) and disk properties. In this paper we have developed 
an iterative 2D, bulge-disk-bar decomposition technique using GALFIT and applied it to 
H-hand images to a complete sample (SI) of 143 bright {Mb < —19.3) low-to-moderately 
inclined {i < 70°) spirals from the OSU Bright Spiral Galaxy Survey (OSUBSGS). The 
sample has primarily spirals with Hubble type SO/a to Sc and stellar mass M^ > 1.0 x IO^^Mq. 
We performed two-component bulge-disk decomposition, as well as three-component bulge- 
disk-bar decomposition on the 2D light distribution of all galaxies, taking into account the 
PSF. We use an exponential profile for the disk, and Sersic profiles for the bulge and bar. A 
number of quantitative indicators, including bar classification from ellipse fits, are used to 
pick either the bulge-disk-bar decomposition or bulge-disk decomposition, as the best final 
fit for a galaxy. Our main results are the following. 

1. We find that it is necessary to include the bar component in 2D decomposition of 
barred galaxies, otherwise, the bulge-to-total ratio (B/T) will be overestimated and 
the disk properties may be skewed. Examples of the effect of including the bar are 
shown for the prominently barred galaxies NGC 4643 (Figure [5l Table [3]) and NGC 
4548 (Figure El Table H. 

2. We find that out of the 143 low-to-moderately inclined {i < 70°) spirals in our sample, 
75 of 143 or ~ 52% are better fit with a Stage 3 bulge-disk-bar decomposition than 
a Stage 2 bulge-disk decomposition. There are also eight galaxies with pure bar-disk 
fits. The resulting H-hand bar fraction, defined as the fraction of disk galaxies that 
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are barred, is 58.0 ± 4.13% (84 of 143). This fraction is in excellent agreement with 
the i7-band bar fraction of 60% reported by MJ07, based on ellipse fits of the same 
OSUBSGS sample, with a more conservative inclination cut {i < 60°). 

//-band images tend to trace the overall mass fairly well and are not overly impacted 
by extinction and age gradients. We therefore assume a constant mass-to-light (M/L) 
in the H-hand for the bulge, disk, and bar, and assume their H-hand light fraction 
is a measure of their mass fraction. For our sample SI of 143 bright {Mb < —19.3) 
low-to-moderately inclined (z < 70°) spirals with a mass-weighted mean Hubble type 
of Sab-Sb, we find that 71.6% of the stellar mass is in disks, 9.6% is in stellar bars and 
18.9% is in bulges (with 10.4% in n > 2 bulges and 8.4% in n < 2 bulges). 

If disks and bars are much younger (e.g., ~ 3 Gyr old ) than bulges (e.g., ~ 12 Gyr 
old), then our prescription would underestimate the true B/T by a factor of < 2. 
On the other hand, if the bulge is younger than the disk and happens to harbor a 
significant young population of massive stars, then our prescription will overestimate 
the true B/T mass ratio, and make our current results on the high fraction of low B/T 
bulges (see point 4 below) even stronger. 



4. We explore the relationship between B/T, bulge Sersic index, and Hubble types (Fig 
& Fig. [T^ . Only a small fraction (~ 1%) of bright spirals have high n > A bulges; such 
bulges lie primarily in SO/a to Sab, and have a large B/T > 0.2. A moderate fraction 
(~ 22%) have intermediate 2 < n < 4 bulges; these exist in barred and unbarred SO/a 
to Sd, and their B/T spans a wide range (0.03 to 0.5) with a mean of 0.29. Finally, a 
strikingly large fraction (^ 7Q%) of bright spirals have low n < 2 bulges; such bulges 
exist in barred and unbarred galaxies across a wide range of Hubble types, and their 
B/T varies from 0.01 to 0.4, with most having B/T < 0.2. 

5. Bulges with very high B/T (> 0.4) exist primarily in galaxies with high mass (M^ > 
6 X IO^^Mq) and early types (SO/a to Sab). Bulges with very low B/T (< 0.1) lie 
primarily in lower mass galaxies with later morphologies (Sb to Sc). As many as 
~ 69% of bright spirals have bulges with B/T < O.2.- these bulges are pervasive and 
exist across the whole spectrum of SO/a to Sd (Figure [H 



6. Modeling bars with 2D decomposition allows us to measure bar properties and the 
bar-to-total ratio (Bar/T), which is a measure of bar strength. There is a wide range 
(~ 0.03 to ~ 0.47) in the individual Bar/T at a given Hubble type. The mean Bar/T 
remains fairly constant with Hubble type from Sa to Sb, but shows a possible weak 
decline by about 0.1 from Sb to Sc (See Figure [19] and Figure [201) . The bar fraction 
(Table [lO]) declines with B/T; it is highest (~ 68%) for bright spirals with B/T < 0.2, 
and lower (~ 36%) by nearly a factor of two in spirals with B/T > 0.2. 
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It may be tempting to infer this to mean the formation of two-thirds of low-B/T bulges 
is related to spontaneous or/and tidally induced bars. Such a cause-effect relationship 
is not the only scenario consistent with this result. It may also be possible that bar 
instabilities are favored in galaxies with low B/T and no inner Lindblad resonances. 
Under these conditions, the swing amplifier model with a feedback loop may be re- 
sponsible for bar formation and partly account for the high bar fraction in galaxies of 
low B/T. 

7. We compare the observed distribution of bulge B/T and n in high mass (M^ > 1.0 x 
IO^^Mq) spirals with predictions from a set of ACDM cosmological semi-analytical 
models (Table [TTl Table [T21 and Figs. [21]to[2ZD. Major mergers are considered as 
those with stellar mass ratio M1/M2 > 1/4. In the models, a high mass spiral can 
have a bulge with a present-day low B/T < 0.2 only if it did not undergo a major 
merger since z < 2 (Fig. [23]) . The model merger history shows that only ~ 13.5% 
and ~ 20% of the high mass spirals experience major mergers since z < 2 and z < 4, 
respectively. The fraction (~ 1.6%) of high mass spirals which have undergone a major 
merger since z < 4 and host a bulge with a present-day low B/T < 0.2 is a factor 
of over thirty smaller than the observed fraction (~ 66%) of high mass spirals with 
B/T < 0.2(Table [TT|) . Thus, bulges built via major mergers since z < 4, over the last 
12 Gyr, seriously fail to account for most of the low B/T < 0.2 bulges present in two- 
thirds of high-mass spirals. The overall picture that emerges is that the observed large 
frequency (~ 66%) of high mass (M^ > 1.0 x 10^°Mq) spirals with low present-day 
B/T < 0.2 can be accounted for in our hierarchical models by high mass spirals, which 
have not undergone a major merger since z < 2, and most of which have not even 
experienced a major merger since z < 4. Most of these present-day low B/T < 0.2 
bulges are likely to have been built by a combination of minor mergers and/or secular 
processes since z < 4. 
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Table 1. Sample SI of Bright {Mb < -19.3) Low-to-Moderately Inclined {i < 70' 

in OSUBSGS (N=143) 
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(Mpc) 


(mag) 


(mag) 


(Me) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


ESO138-10 


PSF+Bulge+Disk 


SA(s)cd 


A 


15.75 


-20.2 


- 


- 


IC0239 


PSF+Bulge+Disk+Bar 


SAB(rs)cd 


AB 


15.21 


-19.3 


0.70 


1.34e+10 


IC4444 


Bulge+Disk 


SAB(rs)bc 


AB 


28.82 


-20.9 


0.64 


4.35e+10 


IC5325 


PSF+Bulge+Disk 


SAB(rs)bc 


AB 


19.39 


-19.6 


0.56 


9.78e+09 


NGC0150 


PSF+Bulge+Disk+Bar 


SB{rs)bc 


B 


20.57 


-20.1 


0.64 


2.17e+10 


NGC0157 


PSF+Bulge+Disk 


SAB(rs)bc 


AB 


22.39 


-21.4 


0.59 


5.81e+10 


NGC0210 


Bulge+Disk+Bar 


SAB(s)b 


AB 


21.75 


-20.4 


0.71 


3.77e+10 


NGC0278 


PSF+Bulge+Disk 


SAB(rs)b 


AB 


12.64 


-19.5 


0.64 


1.24e+10 


NGC0289 


PSF+Bulge+Disk+Bar 


SAB(rs)bc 


AB 


20.79 


-20.2 


0.73 


3.33e+10 


NGC0428 


PSF+Bulge+Disk+Bar 


SAB(s)m 


AB 


15.96 


-19.4 


0.44 


5.28e+09 


NGC0488 


PSF+Bulge+Disk 


SA(r)b 


A 


31.39 


-21.8 


0.87 


2.44e+ll 


NGC0578 


PSF+Bulge+Disk+Bar 


SAB(rs)c 


AB 


20.89 


-20.5 


0.51 


1.87e+10 


NGC0613 


PSF+Bulge+Disk+Bar 


SB(rs)bc 


B 


18.75 


-20.8 


0.68 


4.92e+10 


NGC0685 


Bar+Disk 


SAB(r)c 


AB 


16.29 


-19.4 


0.46 


5.47e+09 


NGC0779 


Bulge+Disk+Bar 


SAB(r)b 


AB 


18.54 


-20.0 


0.79 


3.64C+10 


NGC0864 


PSF+Bulge+Disk+Bar 


SAB(rs)c 


AB 


21.43 


-20.6 


0.55 


2.31e+10 


NGC0908 


PSF+Bulge+Disk 


SA(s)c 


A 


19.07 


-21.3 


0.65 


6.60e+10 


NGC1042 


PSF+Bulge+Disk 


SAB{rs)cd 


AB 


17.89 


-20.2 


0.54 


1.63e+10 


NGC1073 


PSF+Bulge+Disk+Bar 


SB(rs)c 


B 


16.29 


-19.9 


0.50 


1.04e+10 


NGC1084 


PSF+Bulge+Disk 


SA(s)c 


A 


18.32 


-20.6 


0.58 


2.64e+10 


NGC1087 


Bar+Disk 


SAB(rs)c 


AB 


20.36 


-20.6 


0.52 


2.09e+10 


NGC1187 


PSF+Bulge+Disk+Bar 


SB(r)c 


B 


17.46 


-20.2 


0.56 


1.66C+10 


NGC1241 


PSF+Bulge+Disk+Bar 


SB(rs)b 


B 


28.50 


-21.7 


0.85 


2.05e+ll 


NGC1300 


PSF+Bulge+Disk+Bar 


SB(rs)bc 


B 


20.14 


-20.9 


0.68 


5.39e+10 


NGC1302 


PSF+Bulge+Disk+Bar 


{R)SB(r)0 


B 


21.43 


-20.2 


0.89 


6.37e+10 


NGC1309 


PSF+Bulge+Disk 


SA(s)bc 


A 


27.86 


-20.5 


0.44 


1.46e+10 


NGC1317 


Bulge+Disk+Bar 


SAB(r)a 


AB 


18.11 


-20.2 


0.89 


6.42C+10 


NGC1350 


PSF+Bulge+Disk+Bar 


(R')SB(r)ab 


B 


18.11 


-21.1 


0.87 


1.38e+ll 


NGC1371 


PSF+Bulge+Disk+Bar 


SAB(rs)a 


AB 


18.32 


-20.0 


0.90 


5.43e+10 


NGC1385 


PSF+Bar+Disk 


SB(s)cd 


B 


18.75 


-20.3 


0.51 


1.61e+10 


NGC1511 


Bulge+Disk 


SAa;pec 


A 


16.18 


-19.7 


0.57 


1.15e+10 


NGC1559 


PSF+Bar+Disk 


SB(s)cd 


B 


15.32 


-20.4 


0.35 


9.33e+09 


NGC1703 


PSF+Bulge+Disk 


SB(r)b 


B 


18.64 


-19.4 


0.56 


8.43e+09 


NGC1792 


PSF+Bulge+Disk 


SA{rs)bc 


A 


14.57 


-20.5 


0.68 


3.71C+10 


NGC1808 


PSF+Bulge+Disk+Bar 


(R)SAB{s)a 


AB 


11.57 


-19.9 


0.81 


3.57e+10 


NGC1964 


PSF+Bulge+Disk+Bar 


SAB(s)b 


AB 


21.43 


-20.8 


0.77 


6.54e+10 


NGC2090 


Bulge+Disk+Bar 


SA{rs)c 


A 


10.93 


-19.3 


0.79 


1.90C+10 


NGC2139 


Bulge+Disk 


SAB{rs)cd 


AB 


24.00 


-20.3 


0.36 


8.35e+09 


NGC2196 


PSF+Bulge+Disk 


(R')SA(s)a 


A 


30.86 


-21.0 


0.81 


9.22e+10 


NGC2442 


PSF+Bulge+Disk+Bar 


SAB(s)bc;pcc 


AB 


18.32 


-20.8 


0.82 


8.28e+10 


NGC2559 


Bulge+Disk+Bar 


SB(s)bc;pec 


B 


21.43 


-21.0 


- 


- 


NGC2566 


Bulge+Disk 


(R')SB{rs)ab;pec 


B 


22.61 


-20.7 


0.81 


6.98e+10 


NGC2775 


PSF+Bulge+Disk 


SA(r)ab 


A 


18.21 


-20.6 


0.90 


9.46e+10 


NGC3059 


Bar+Disk 


SB(rs)c 


B 


15.86 


-20.1 


0.68 


2.52e+10 
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Table 1 — Continued 



Galaxy Name 


Best Fit 


Hubble Type 


Bar Type 


D 


Mb 


B~V 


M* 






(RC3) 


(RC3) 


(Mpc) 


(mag) 


(mag) 


{Mq) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


NGC3166 


Bulge+Disk+Bar 


SAB(rs)0 


AB 


23.57 


-20.2 


0.93 


7.41e+10 


NGC3169 


PSF+Bulge+Disk 


SA(s)a;pec 


A 


21.11 


-20.3 


0.85 


6.09e+10 


NGC3223 


PSF+Bulge+Disk 


SA(s)b 


A 


40.82 


-21.9 


0.82 


2.36e+ll 


NGC3227 


PSF+Bulge+Disk+Bar 


SAB(s)a;pec 


AB 


22.07 


-20.1 


0.82 


4.42e+10 


NGC3261 


PSF+Bulgc+Disk+Bar 


SB(rs)b 


B 


35.79 


-21.3 


- 


- 


NGC3275 


PSF+Bulgc+Disk+Bar 


SB(r)ab 


B 


45.43 


-21.2 


- 


- 


NGC3319 


PSF+Bar+Disk 


SB(rs)cd 


B 


12.32 


-19.4 


0.41 


4.76e+09 


NGC3338 


PSF+Bulge+Disk+Bar 


SA(s)c 


A 


24.43 


-20.6 


0.59 


2.67e+10 


NGC3423 


Bulge+Disk 


SA(s)cd 


A 


11.68 


-19.5 


0.45 


5.89e+09 


NGC3504 


Bulge+Disk+Bar 


(R)SAB(s)ab 


AB 


21.43 


-20.5 


0.72 


4.06e+10 


NGC3513 


Bulge+Disk+Bar 


SB(rs)c 


B 


18.21 


-19.9 


0.43 


7.52C+09 


NGC3583 


Bulge+Disk+Bar 


SB(s)b 


B 


36.43 


-20.8 


- 


- 


NGC3596 


PSF+Bulge+Disk 


SAB(rs)c 


AB 


24.64 


-19.6 


- 


- 


NGC3646 


PSF+Bulge+Disk 


Ring 


- 


59.79 


-22.9 


0.65 


2.82e+ll 


NGC3675 


PSF+Bulge+Disk 


SA(s)b 


A 


13.71 


-20.2 


- 


- 


NGC3684 


PSF+Bulge+Disk 


SA(rs)bc 


A 


25.07 


-19.4 


0.62 


l.OOe+10 


NGC3686 


PSF+Bulge+Disk+Bar 


SB(s)bc 


B 


25.18 


-19.6 


0.57 


9.97e+09 


NGC3705 


PSF+Bulge+Disk+Bar 


SAB(r)ab 


AB 


18.21 


-19.9 


0.79 


3.15e+10 


NGC3726 


PSF+Bulge+Disk+Bar 


SAB(r)c 


AB 


18.21 


-20.6 


0.49 


1.93e+10 


NGC3810 


PSF+Bulge+Disk 


SA{rs)c 


A 


18.11 


-20.1 


0.58 


1.72e+10 


NGC3885 


PSF+Bulge+Disk 


SA{s)0 


A 


29.79 


-19.6 


0.95 


4.62e+10 


NGC3887 


PSF+Bulge+Disk+Bar 


SB(r)bc 


B 


20.68 


-20.4 


- 


- 


NGC3893 


PSF+Bulge+Disk 


SAB(rs)c 


AB 


18.21 


-19.9 


- 


- 


NGC3938 


PSF+Bulge+Disk 


SA(s)c 


A 


18.21 


-20.0 


0.52 


1.23e+10 


NGC3949 


PSF+Bulge+Disk 


SA(s)bc 


AB 


18.21 


-19.9 


0.45 


8.66e+09 


NGC4027 


Bar+Disk 


SB(s)dm 


B 


27.43 


-20.6 


0.54 


2.25e+10 


NGC4030 


PSF+Bulge+Disk 


SA(s)bc 


A 


27.75 


-20.8 


- 


- 


NGC4051 


PSF+Bulge+Disk+Bar 


SAB(rs)bc 


AB 


18.21 


-20.0 


0.65 


1.95e+10 


NGC4062 


PSF+Bulge+Disk 


SA(s)c 


A 


10.39 


-19.5 


0.76 


2.07e+10 


NGC4123 


PSF+Bulge+Disk+Bar 


SB(r)c 


B 


17.68 


-19.8 


0.61 


1.50e+10 


NGC4145 


PSF+Bulge+Disk 


SAB(rs)d 


AB 


22.18 


-20.1 


0.51 


1.34e+10 


NGC4151 


PSF+Bulge+Disk+Bar 


(R')SAB(rs)ab 


AB 


21.75 


-20.1 


0.73 


2.93e+10 


NGC4212 


PSF+Bulge+Disk 


SAc 


A 


18.00 


-20.4 


0.67 


3.28e+10 


NGC4254 


PSF+Bulge+Disk 


SA{s)c 


A 


18.00 


-22.6 


0.57 


1.61e+ll 


NGC4293 


PSF+Bulge+Disk+Bar 


(R)SB(s)0 


B 


18.21 


-20.1 


0.90 


5.94e+10 


NGC4303 


Bulge+Disk+Bar 


SAB(rs)bc 


AB 


16.29 


-21.8 


0.53 


6.76e+10 


NGC4314 


PSF+Bulge+Disk+Bar 


SB(rs)a 


B 


10.39 


-19.8 


0.85 


3.69e+10 


NGC4394 


PSF+Bulge+Disk+Bar 


(R)SB(r)b 


B 


18.00 


-19.4 


0.85 


2.61e+10 


NGC4414 


PSF+Bulge+Disk+Bar 


SA(rs)c 


A 


10.39 


-20.0 


0.84 


4.38e+10 


NGC4450 


PSF+Bulge+Disk+Bar 


SA(s)ab 


A 


18.00 


-21.9 


0.82 


2.22e+ll 


NGC4487 


PSF+Bulge+Disk+Bar 


SAB(rs)cd 


AB 


21.32 


-19.6 


- 


- 


NGC4490 


PSF+Bulge+Disk 


SB(s)d;pec 


B 


8.36 


-21.7 


0.43 


4.10e+10 


NGC4527 


Bulge+Disk+Bar 


SAB(s)bc 


AB 


14.46 


-21.5 


0.86 


1.87C+11 


NGC4548 


PSF+Bulge+Disk+Bar 


SB(rs)b 


B 


18.00 


-20.8 


0.81 


7.85e+10 
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Table 1 — Continued 



Galaxy Name 


Best Fit 


Hubble Type 


Bar Type 


D 


Mb 


B~V 


M* 






(RC3) 


(RC3) 


(Mpc) 


(mag) 


(mag) 


(Mq) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


NGC4593 


PSF+Bulge+Disk+Bar 


(R)SB{rs)b 


B 


42.32 


-20.8 


_ 


_ 


NGC4618 


PSF+Bar+Disk 


SB(rs)m 


B 


7.82 


-19.3 


0.44 


4.62e+09 


NGC4643 


Bulge+Disk+Bar 


SB(rs)0 


B 


27.54 


-19.9 


0.96 


6.49e+10 


NGC4647 


Bulge+Disk 


SAB(rs)c 


AB 


18.00 


-19.8 


0.65 


1.67e+10 


NGC4651 


PSF+Bulge+Disk 


SA(rs)c 


A 


18.00 


-19.6 


0.57 


1.05e+10 


NGC4654 


Bulge+Disk 


SAB(rs)cd 


AB 


18.00 


-20.6 


0.60 


2.99e+10 


NGC4665 


PSF+Bulge+Disk+Bar 


SB(s)0 


B 


19.18 


-19.3 


- 


- 


NGC4689 


PSF+Bulge+Disk 


SA(rs)bc 


A 


18.00 


-20.7 


0.65 


3.84e+10 


NGC4691 


PSF+Bulge+Disk+Bar 


{R)SB(s)0;pec 


B 


24.11 


-19.6 


0.58 


1.04e+10 


NGC4698 


PSF+Bulge+Disk 


SA(s)ab 


A 


18.00 


-19.9 


0.91 


5.15e+10 


NGC4699 


Bulge+Disk+Bar 


SAB(rs)b 


AB 


27.54 


-21.5 


0.89 


2.06e+ll 


NGC4772 


PSF+Bulge+Disk+Bar 


SA(s)a 


A 


17.46 


-19.4 


0.92 


3.30e+10 


NGC4775 


PSF+Bulge+Disk 


SA(s)d 


A 


28.50 


-20.4 


- 


- 


NGC4781 


PSF+Bulge+Disk 


SB(rs)d 


B 


24.11 


-20.8 


- 


- 


NGC4818 


PSF+Bulge+Disk+Bar 


SAB(rs)ab;pcc 


AB 


23.04 


-19.7 


0.89 


3.88e+10 


NGC4856 


PSF+Bulge+Disk+Bar 


SB(s)0 


B 


22.61 


-20.2 


0.99 


9.04e+10 


NGC4902 


PSF+Bulge+Disk+Bar 


SB(r)b 


B 


42.00 


-21.4 


0.69 


8.33e+10 


NGC4930 


PSF+Bulge+Disk+Bar 


SB(rs)b 


B 


37.50 


-21.2 


0.90 


1.61e+ll 


NGC4939 


PSF+Bulge+Disk 


SA{s)bc 


A 


47.46 


-22.2 


0.64 


1.43e+ll 


NGC4941 


PSF+Bulge+Disk 


(R)SAB(r)ab 


AB 


6.86 


-19.4 


0.84 


2.50e+10 


NGC4995 


PSF+Bulge+Disk+Bar 


SAB(rs)b 


AB 


30.00 


-20.6 


0.87 


8.25e+10 


NGC5054 


Bulge+Disk 


SA(s)bc 


A 


29.25 


-21.1 


0.76 


8.31e+10 


NGC5085 


Bulge+Disk 


SA{s)c 


A 


30.96 


-19.3 


0.87 


2.54e+10 


NGC5101 


PSF+Bulge+Disk+Bar 


(R)SB(rs)0 


B 


29.36 


-20.8 


1.00 


1.69e+ll 


NGC5121 


PSF+Bulge+Disk 


(R')SA(s)a 


A 


23.68 


-19.4 


0.95 


3.68e+10 


NGC5161 


PSF+Bulge+Disk 


SA(s)c 


A 


35.89 


-21.7 


0.79 


1.65e+ll 


NGC5247 


PSF+Bulge+Disk 


SA(s)bc 


A 


23.79 


-21.2 


0.54 


3.86e+10 


NGC5371 


PSF+Bulge+Disk+Bar 


SAB(rs)bc 


AB 


40.50 


-22.1 


0.70 


1.74e+ll 


NGC5427 


PSF+Bulge+Disk 


SA(s)c;pec 


A 


40.82 


-21.2 


0.57 


4.61e+10 


NGC5483 


Bulge+Disk+Bar 


SA(s)c 


A 


26.46 


-20.3 


- 


- 


NGC5643 


PSF+Bulge+Disk+Bar 


SAB(rs)c 


AB 


18.11 


-20.9 


0.74 


6.68e+10 


NGC5676 


PSF+Bulge+Disk 


SA(rs)bc 


A 


36.96 


-21.5 


0.68 


9.01e+10 


NGC5701 


PSF+Bulge+Disk+Bar 


{R)SB(rs)0 


B 


27.96 


-19.9 


0.88 


4.67e+10 


NGC5713 


PSF+Bulge+Disk 


SAB(rs)bc;pec 


AB 


32.57 


-21.1 


0.64 


5.21e+10 


NGC5850 


PSF+Bulge+Disk+Bar 


SB(r)b 


B 


30.54 


-21.5 


0.79 


1.37e+ll 


NGC5921 


PSF+Bulge+Disk+Bar 


SB(r)bc 


B 


27.00 


-20.6 


0.66 


3.51e+10 


NGC5962 


Bulge+Disk+Bar 


SA(r)c 


A 


34.07 


-20.9 


0.64 


4.43e+10 


NGC6215 


PSF+Bulge+Disk 


SA(s)c 


A 


21.96 


-20.4 


0.54 


1.84e+10 


NGC6221 


Bulge+Disk+Bar 


SB(s)bc;pec 


B 


20.79 


-21.6 


0.74 


1.31e+ll 


NGC6300 


PSF+Bulge+Disk+Bar 


SB(rs)b 


B 


15.32 


-20.3 


0.78 


4.49e+10 


NGC6384 


Bulge+Disk 


SAB(r)bc 


AB 


28.50 


-21.5 


0.72 


1.05e+ll 


NGC6753 


PSF+Bulge+Disk+Bar 


(R)SA(r)b 


A 


43.82 


-21.6 


0.83 


1.81e+ll 


NGC6782 


PSF+Bulge+Disk+Bar 


{R)SAB{r)a 


AB 


54.43 


-21.4 


- 


- 


NGC6902 


PSF+Bulge+Disk 


SA(r)b 


A 


38.25 


-21.2 


0.71 


8.12e+10 
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Table 1 — Continued 



Galaxy Name 


Best Fit 


Hubble Type 


Bar Type 


D 


Mb 


B-V 


M* 






(RCS) 


(RC3) 


(Mpc) 


(mag) 


(mag) 


(Mq) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


NGC6907 


PSF+Bulgc+Disk+Bar 


SB(s)bc 


B 


46.07 


-21.8 


0.69 


1.21e+ll 


NGC7083 


PSF+Bulge+Disk 


SA(s)bc 


A 


41.46 


-21.8 


0.65 


1.05e+ll 


NGC7205 


PSF+Bulge+Disk 


SA{s)bc 


A 


21.96 


-20.7 


0.60 


3.11e+10 


NGC7213 


PSF+Bulge+Disk 


SA{s)a 


A 


23.57 


-21.0 


0.89 


1.26e+ll 


NGC7217 


PSF+Bulge+Disk 


(R)SA(r)ab 


A 


17.14 


-20.5 


0.90 


8.38e+10 


NGC7412 


PSF+Bulge+Disk+Bar 


SB(s)b 


B 


22.61 


-20.2 


0.53 


1.50C+10 


NGC7479 


PSF+Bulge+Disk+Bar 


SB(s)c 


B 


34.71 


-21.7 


0.75 


1.44e+ll 


NGC7552 


Bulge+Disk+Bar 


(R')SB{s)ab 


B 


20.89 


-20.5 


0.68 


3.49e+10 


NGC7723 


PSF+Bulge+Disk+Bar 


SB(r)b 


B 


25.39 


-20.5 


0.73 


4.32e+10 


NGC7727 


PSF+Bulge+Disk 


SAB(s)a;pec 


AB 


24.96 


-20.8 


0.91 


1.13e+ll 


NGC7741 


Bulge+Disk+Bar 


SB(s)cd 


B 


13.18 


-19.3 


0.53 


6.67e+09 



Note. — Columns are: (1) Galaxy name. (2) The best fit chosen based on the criteria outlined in § 13.31 (3) 
Hubble type from RC3 (de Vaucouleurs et al. 1991). (4) RC3 bar type, which is based on visual inspection of 
optical images and runs as 'B'='strongly barred', 'AB'='weakly barred', and 'A'='unbarred'. (5) The distance in 
Mpc from the Nearby Galaxies Catalog (TuUy 1988) scaled to a Hubble constant of 70 km s~^ Mpc~^. Distances 
are adjusted assuming the Galaxy is retarded by 300 km s~^ from the universal expansion by the mass of the Virgo 
Cluster. Exceptions are NGC 6753, NGC 6782, and NGC 6907 where distances are from RC3, NGC 3504 where 
the distance is from Kenney et al. (1993), and NGC 4314 where the distance is from Benedict et al. (1996). (6) 
Absolute B-band magnitude from Hyperleda. Corrections have been applied for galactic and internal extinction. 
Galactic absorption is computed based on Schlegel et al. (1998). Internal extinction is corrected for following 
Bottinelli et al. (1995). K-correction co-efficients are drawn from RC2 (de Vaucouleurs et al. 1976). (7) B — V 
color from Hyperleda. Corrections for extinction and k-correction have been applied as in column 6. (8) Stellar 
mass, calculated as outlined in ^ 12.21 



Table 2. Parameters From 2D Luminosity Decomposition for Sample SI (N=143) 



Galaxy Name 


Best Fit 


B/T 


D/T 


Bar/T 


B/D 


Bulge re 


Bulge re 


Bulge n 


Disk h 


Diskh 


Bar re 


Bar re 


Bar n 






% 


% 


% 




// 


kpc 




// 


kpc 


// 


kpc 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


ESO138-10 


PSF+Bulgc+Disk 


16.25 


83.75 


_ 


0.19 


10.94 


0.86 


2.07 


27.10 


2.14 


_ 


_ 


_ 


IC0239 


PSF+Bulgc+Disk+Bar 


2.45 


92.65 


4.91 


0.03 


6.63 


0.41 


0.23 


34.56 


2.15 


13.95 


0.87 


0.19 


IC4444 


Bulgc+Disk 


31.97 


68.03 


- 


0.47 


7.94 


1.06 


2.33 


17.67 


2.36 


- 


- 


- 


IC5325 


PSF+Bulgc+Disk 


6.79 


93.21 


- 


0.07 


15.17 


1.56 


1.75 


23.18 


2.39 


- 


- 


- 


NGC0150 


PSF+Bulge+Disk+Bar 


5.72 


81.18 


13.11 


0.07 


4.22 


0.44 


0.05 


33.57 


3.52 


17.82 


1.87 


0.39 


NGC0157 


PSF+Bulge+Disk 


2.29 


97.71 


- 


0.02 


3.44 


0.41 


0.44 


31.82 


3.78 


- 


- 


- 


NGC0210 


Bulge+Disk+Bar 


30.07 


48.54 


21.39 


0.62 


4.81 


0.57 


1.71 


83.21 


9.92 


29.41 


3.51 


0.31 


NGC0278 


PSF+Bulge+Disk 


4.23 


95.77 


- 


0.04 


3.46 


0.15 


0.57 


13.48 


0.58 


- 


- 


- 


NGC0289 


PSF+Bulge+Disk+Bar 


8.83 


85.99 


5.18 


0.10 


4.91 


0.57 


0.46 


20.28 


2.36 


18.90 


2.19 


0.04 


NGC0428 


PSF+Bulge+Disk+Bar 


5.50 


71.64 


22.86 


0.08 


9.27 


0.67 


0.38 


32.55 


2.34 


30.70 


2.21 


0.37 


NGC0488 


PSF+Bulge+Disk 


21.72 


78.28 


- 


0.28 


9.91 


1.52 


3.07 


38.94 


5.96 


- 


- 


- 


NGC0578 


PSF+Bulge+Disk+Bar 


2.05 


93.63 


4.32 


0.02 


4.19 


0.46 


0.56 


40.73 


4.47 


14.95 


1.64 


0.18 


NGC0613 


PSF+Bulge+Disk+Bar 


13.14 


56.97 


29.90 


0.23 


6.12 


0.64 


1.05 


44.69 


4.64 


62.49 


6.49 


0.53 


NGC0685 


Bar+Disk 


- 


96.53 


3.47 


- 


- 


- 


- 


41.24 


4.28 


20.68 


2.15 


0.16 


NGC0779 


Bulge+Disk+Bar 


16.18 


57.69 


26.13 


0.28 


6.04 


0.59 


2.31 


38.74 


3.79 


29.00 


2.84 


0.25 


NGC0864 


PSF+Bulge+Disk+Bar 


2.77 


86.88 


10.35 


0.03 


3.52 


0.38 


0.66 


39.54 


4.22 


21.07 


2.25 


0.37 


NGC0908 


PSF+Bulge+Disk 


8.86 


91.14 


- 


0.10 


7.18 


0.84 


1.62 


48.96 


5.72 


- 


- 


- 


NGC1042 


PSF+Bulge+Disk 


3.00 


97.00 


- 


0.03 


5.91 


0.57 


0.09 


39.83 


3.86 


- 


- 


- 


NGC1073 


PSF+Bulge+Disk+Bar 


2.36 


79.58 


18.06 


0.03 


5.47 


0.46 


0.41 


50.02 


4.17 


34.24 


2.85 


0.83 


NGC1084 


PSF+Bulge+Disk 


5.23 


94.77 


- 


0.06 


3.24 


0.32 


0.77 


18.68 


1.82 


- 


- 


- 


NGC1087 


Bar+Disk 


- 


93.13 


6.87 


- 


- 


- 


- 


28.55 


2.78 


7.79 


0.76 


1.27 


NGC1187 


PSF+Bulge+Disk+Bar 


6.94 


64.96 


28.10 


0.11 


4.81 


0.51 


1.57 


38.20 


4.06 


51.85 


5.51 


0.72 


NGC1241 


PSF+Bulge+Disk+Bar 


11.70 


69.10 


19.21 


0.17 


2.61 


0.70 


0.91 


27.69 


7.43 


17.89 


4.80 


0.49 


NGC1300 


PSF+Bulge+Disk+Bar 


11.04 


74.94 


14.02 


0.15 


4.63 


0.51 


1.19 


63.53 


6.96 


65.89 


7.21 


0.25 


NGC1302 


PSF+Bulge+Disk+Bar 


17.65 


63.79 


18.57 


0.28 


4.96 


0.59 


2.01 


44.92 


5.34 


20.44 


2.43 


0.76 


NGC1309 


PSF+Bulge+Disk 


27.07 


72.93 


- 


0.37 


10.86 


1.68 


2.27 


16.67 


2.58 


- 


- 


- 


NGC1317 


Bulge+Disk+Bar 


13.86 


41.48 


44.65 


0.33 


3.99 


0.53 


1.94 


47.37 


6.31 


20.22 


2.69 


2.19 


NGC1350 


PSF+Bulge+Disk+Bar 


19.52 


69.99 


10.50 


0.28 


7.41 


0.95 


1.59 


52.77 


6.73 


51.07 


6.51 


0.31 


NGC1371 


PSF+Bulge+Disk+Bar 


12.46 


80.16 


7.38 


0.16 


5.17 


0.52 


1.11 


29.34 


2.93 


18.92 


1.89 


0.36 


NGC1385 


PSF+Bar+Disk 


- 


75.54 


24.46 


- 


- 


- 


- 


28.69 


2.87 


15.63 


1.56 


1.28 


NGC1511 


Bulge+Disk 


28.47 


71.53 


- 


0.40 


27.78 


2.55 


1.35 


23.99 


2.20 


- 


- 


- 


NGC1559 


PSF+Bar+Disk 


- 


97.06 


2.94 


- 


- 


- 


- 


31.16 


2.86 


9.84 


0.90 


0.53 


NGC1703 


PSF+Bulge+Disk 


6.32 


93.68 


- 


0.07 


2.83 


0.30 


0.56 


18.24 


1.91 


- 


- 


- 
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Galaxy Name 


Best Fit 


B/T 


D/T 


Bar/T 


B/D 


Bulge re 


Bulge re 


Bulge n 


Disk h 


Diskh 


Bar Te. 


Bar re 


Bar n 






% 


% 


% 




It 


kpc 




II 


kpc 


II 


kpc 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 



NGC1792 PSF+Bulge+Disk 2.65 97.35 

NGC1808 PSF+Bulge+Disk+Bar 3.92 74.27 

NGC1964 PSF+Bulge+Disk+Bar 41.78 49.87 

NGC2090 Bulge+Disk+Bar 14.67 58.39 

NGC2139 Bulge+Disk 15.16 84.84 

NGC2196 PSF+Bulge+Disk 46.38 53.62 

NGC2442 PSF+Bulge+Disk+Bar 13.33 62.47 

NGC2559 Bulge+Disk+Bar 6.98 81.64 

NGC2566 Bulge+Disk 22.71 77.29 

NGC2775 PSF+Bulge+Disk 60.87 39.13 

NGC3059 Bar+Disk - 90.68 

NGC3166 Bulge+Disk+Bar 24.97 50.74 

NGC3169 PSF+Bulge+Disk 38.53 61.47 

NGC3223 PSF+Bulge+Disk 11.71 88.29 

NGC3227 PSF+Bulge+Disk+Bar 10.78 45.23 

NGC3261 PSF+Bulge+Disk+Bar 15.06 70.15 

NGC3275 PSF+Bulge+Disk+Bar 14.15 63.44 

NGC3319 PSF+Bar+Disk - 91.41 

NGC3338 PSF+Bulge+Disk+Bar 5.29 81.96 

NGC3423 Bulge+Disk 28.10 71.90 

NGC3504 Bulge+Disk+Bar 23.15 39.49 

NGC3513 Bulge+Disk+Bar 4.45 92.08 

NGC3583 Bulge+Disk+Bar 8.62 53.40 

NGC3596 PSF+Bulge+Disk 12.15 87.85 

NGC3646 PSF+Bulge+Disk 13.53 86.47 

NGC3675 PSF+Bulge+Disk 44.42 55.58 

NGC3684 PSF+Bulge+Disk 26.92 73.08 

NGC3686 PSF+Bulge+Disk+Bar 3.02 95.00 

NGC3705 PSF+Bulge+Disk+Bar 27.73 59.20 

NGC3726 PSF+Bulge+Disk+Bar 3.09 95.10 

NGC3810 PSF+Bulge+Disk 38.43 61.57 

NGC3885 PSF+Bulge+Disk 27.10 72.90 

NGC3887 PSF+Bulge+Disk+Bar 4.19 83.01 



21.82 
8.35 
26.94 



24.19 
11.37 



9.32 
24.29 



44.00 
14.79 
22.41 
8.59 
12.75 

37.37 
3.47 
37.98 



1.98 
13.07 
1.81 



12.80 



0.03 
0.05 
0.84 
0.25 
0.18 
0.86 
0.21 
0.09 
0.29 
1.56 

0.49 
0.63 
0.13 
0.24 
0.21 
0.22 

0.06 
0.39 
0.59 
0.05 
0.16 
0.14 
0.16 
0.80 
0.37 
0.03 
0.47 
0.03 
0.62 
0.37 
0.05 



3.72 

2.78 
8.31 
17.68 
8.27 
13.42 
4.78 
4.01 
1.74 
47.35 

3.25 
11.01 
5.39 
4.84 
2.95 
2.21 

4.39 
24.09 
2.37 
10.73 
1.29 
5.34 
3.60 
28.41 
11.43 
4.39 
8.62 
5.64 
12.35 
3.13 
5.27 



0.31 
0.19 
0.97 
1.12 
1.04 
2.12 
0.46 
0.43 
0.20 
4.37 

0.30 
0.96 
1.06 
0.44 
0.52 
0.48 

0.39 
1.39 
0.25 
0.88 
0.19 
0.43 
1.04 
1.42 
1.10 
0.31 
0.63 
0.37 
0.82 
0.41 
0.44 



0.84 
0.45 
2.58 
2.19 
1.53 
2.38 
1.68 
1.38 
4.42 
4.85 

0.81 
2.41 
1.74 
0.32 
1.96 
1.85 

0.76 
2.39 
1.04 
1.46 
1.28 
0.67 
1.20 
2.66 
2.15 
0.66 
1.83 
0.53 
1.28 
0.46 
1.01 



38.50 
38.78 
18.59 
160.02 
18.12 
28.38 
78.46 
34.65 
23.05 
28.09 
64.22 
20.14 
55.12 
29.05 
44.67 
28.89 
28.25 
66.08 
36.42 
31.21 
29.55 
33.18 
21.78 
18.54 
26.01 
47.12 
15.96 
27.57 
37.77 
49.58 
23.78 
11.94 
42.69 



3.23 

2.71 
2.17 
10.17 
2.27 
4.48 
7.55 
3.72 
2.60 
2.59 
5.92 
1.84 
4.79 
5.73 
4.07 
5.08 
6.20 
14.49 
3.25 
1.80 
3.09 
2.72 
3.19 
1.50 
7.54 
2.35 
1.53 
1.96 
2.74 
3.24 
1.57 
1.57 
3.55 



5.71 
8.05 
33.74 



60.15 
21.40 



23.94 
13.85 



39.51 
14.19 
20.75 
14.60 
20.62 

26.02 
18.32 
14.89 



21.32 
30.00 
41.09 



37.62 



0.40 
0.94 
2.14 



5.79 
2.30 



2.21 
1.26 



3.60 

2.49 
4.55 
3.20 

1.84 

2.72 
1.50 
2.18 



1.52 
2.18 
2.69 



3.13 



0.67 
0.05 
0.43 



0.24 
0.24 



1.43 
0.53 



1.27 
1.05 
0.94 
0.43 
0.39 

0.78 
0.08 
0.68 



0.03 
0.07 
0.05 



0.62 



00 
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Galaxy Name 


Best Fit 


B/T 


D/T 


Bar/T 


B/D 


Bulge re 


Bulge Te 


Bulge n 


Disk h 


Diskh 


Bar re 


Bar re 


Bar n 








% 


% 


% 




// 


kpc 




II 


kpc 


II 


kpc 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 




NGC3893 


PSF+Bulge+Disk 


54.61 


45.39 


_ 


1.20 


20.05 


1.31 


2.05 


23.40 


1.52 


_ 


_ 


_ 




NGC3938 


PSF+Bulge+Disk 


7.30 


92.70 


- 


0.08 


6.43 


0.34 


1.18 


32.84 


1.75 


- 


- 


- 




NGC3949 


PSF+Bulge+Disk 


7.75 


92.25 


- 


0.08 


4.67 


0.22 


0.64 


15.22 


0.72 


- 


- 


- 




NGC4027 


Bar+Disk 


- 


79.48 


20.52 


- 


- 


- 


- 


43.11 


2.03 


26.50 


1.25 


1.80 




NGC4030 


PSF+Bulge+Disk 


47.29 


52.71 


- 


0.90 


17.18 


1.73 


2.22 


20.67 


2.08 


- 


- 


- 




NGC4051 


PSF+Bulge+Disk+Bar 


14.99 


46.99 


38.02 


0.32 


11.29 


0.54 


1.20 


65.81 


3.12 


69.69 


3.31 


0.40 




NGC4062 


PSF+Bulge+Disk 


1.69 


98.31 


- 


0.02 


3.62 


0.19 


0.50 


34.98 


1.79 


- 


- 


- 




NGC4123 


PSF+Bulge+Disk+Bar 


6.87 


79.23 


13.90 


0.09 


6.55 


0.56 


0.59 


56.76 


4.84 


42.98 


3.66 


0.46 




NGC4145 


PSF+Bulge+Disk 


7.75 


92.25 


- 


0.08 


12.90 


0.82 


0.68 


63.59 


4.05 


- 


- 


- 




NGC4151 


PSF+Bulge+Disk+Bar 


41.32 


50.26 


8.42 


0.82 


9.36 


0.62 


0.44 


37.74 


2.49 


57.00 


3.76 


0.10 




NGC4212 


PSF+Bulge+Disk 


3.67 


96.33 


- 


0.04 


1.44 


0.01 


1.38 


22.56 


0.13 


- 


- 


- 




NGC4254 


PSF+Bulge+Disk 


38.98 


61.02 


- 


0.64 


30.46 


5.03 


2.68 


34.19 


5.64 


- 


- 


- 




NGC4293 


PSF+Bulge+Disk+Bar 


8.49 


61.70 


29.81 


0.14 


8.00 


0.40 


0.99 


58.26 


2.88 


55.17 


2.73 


0.46 


1 


NGC4303 


Bulge+Disk+Bar 


8.28 


82.81 


8.91 


0.10 


2.67 


0.30 


1.55 


44.58 


4.93 


32.35 


3.58 


0.55 


CO 


NGC4314 


PSF+Bulge+Disk+Bar 


26.49 


52.44 


21.07 


0.51 


11.22 


0.74 


2.05 


62.35 


4.14 


48.98 


3.25 


0.35 


NGC4394 


PSF+Bulge+Disk+Bar 


15.20 


69.58 


15.22 


0.22 


4.20 


0.22 


1.63 


41.74 


2.22 


28.70 


1.53 


0.59 




NGC4414 


PSF+Bulge+Disk+Bar 


11.48 


83.86 


4.65 


0.14 


3.50 


0.17 


1.72 


21.27 


1.05 


17.66 


0.87 


0.05 




NGC4450 


PSF+Bulge+Disk+Bar 


16.86 


71.49 


11.65 


0.24 


8.19 


1.10 


2.26 


47.16 


6.33 


36.13 


4.85 


0.33 




NGC4487 


PSF+Bulge+Disk+Bar 


3.10 


93.77 


3.13 


0.03 


6.66 


0.47 


0.83 


27.92 


1.99 


14.52 


1.04 


0.27 




NGC4490 


PSF+Bulge+Disk 


6.64 


93.36 


- 


0.07 


21.41 


0.88 


0.35 


44.02 


1.81 


- 


- 


- 




NGC4527 


Bulge+Disk+Bar 


18.24 


63.84 


17.91 


0.29 


5.17 


0.61 


1.99 


55.52 


6.59 


41.19 


4.89 


0.50 




NGC4548 


PSF+Bulge+Disk+Bar 


12.96 


68.65 


18.39 


0.19 


6.97 


0.23 


1.56 


58.22 


1.96 


44.91 


1.51 


0.51 




NGC4593 


PSF+Bulge+Disk+Bar 


25.13 


28.15 


46.72 


0.89 


5.27 


0.97 


1.04 


33.07 


6.11 


41.94 


7.75 


0.64 




NGC4618 


PSF+Bar+Disk 


- 


86.11 


13.89 


- 


- 


- 


- 


46.22 


8.54 


17.05 


3.15 


0.66 




NGC4643 


Bulge+Disk+Bar 


24.96 


54.10 


20.95 


0.46 


5.43 


0.52 


2.53 


48.22 


4.64 


21.30 


2.05 


0.62 




NGC4647 


Bulge+Disk 


33.39 


66.61 


- 


0.50 


12.90 


1.26 


1.61 


42.63 


4.17 


- 


- 


- 




NGC4651 


PSF+Bulge+Disk 


41.67 


58.33 


- 


0.71 


18.52 


1.01 


0.63 


30.46 


1.66 


- 


- 


- 




NGC4654 


Bulge+Disk 


1.74 


98.26 


- 


0.02 


3.78 


0.27 


0.93 


30.16 


2.18 


- 


- 


- 




NGC4665 


PSF+Bulge+Disk+Bar 


14.74 


67.38 


17.89 


0.22 


6.92 


0.37 


1.12 


50.48 


2.73 


30.68 


1.66 


0.70 




NGC4689 


PSF+Bulge+Disk 


5.22 


94.78 


- 


0.06 


6.97 


0.73 


1.14 


41.85 


4.38 


- 


- 


- 




NGC4691 


PSF+Bulge+Disk+Bar 


15.80 


69.09 


15.11 


0.23 


11.62 


0.89 


0.61 


33.90 


2.59 


24.99 


1.91 


0.42 




NGC4698 


PSF+Bulge+Disk 


21.72 


78.28 


- 


0.28 


6.78 


0.48 


1.78 


28.34 


2.02 


- 


- 


- 




NGC4699 


Bulge+Disk+Bar 


19.93 


76.48 


3.59 


0.26 


2.62 


0.27 


2.08 


15.93 


1.62 


13.49 


1.37 


0.02 
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Galaxy Name 


Best Fit 


B/T 


D/T 


Bar/T 


B/D 


Bulge Ve 


Bulge re 


Bulge n 


Disk h 


Diskh 


Bar Te. 


Bar re 


Bar n 






% 


% 


% 




II 


kpc 




II 


kpc 


II 


kpc 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 



NGC4772 PSF+Bulge+Disk+Bar 34.25 51.55 

NGC4775 PSF+Bulge+Disk 32.85 67.15 

NGC4781 PSF+Bulge+Disk 10.18 89.82 

NGC4818 PSF+Bulge+Disk+Bar 7.04 73.74 

NGC4856 PSF+Bulge+Disk+Bar 27.44 61.22 

NGC4902 PSF+Bulge+Disk+Bar 5.59 84.44 

NGC4930 PSF+Bulge+Disk+Bar 23.74 61.29 

NGC4939 PSF+Bulge+Disk 21.37 78.63 

NGC4941 PSF+Bulge+Disk 15.08 84.92 

NGC4995 PSF+Bulge+Disk+Bar 6.62 87.18 

NGC5054 Bulge+Disk 9.64 90.36 

NGC5085 Bulge+Disk 6.04 93.96 

NGC5101 PSF+Bulge+Disk+Bar 32.63 48.31 

NGC5121 PSF+Bulge+Disk 36.43 63.57 

NGC5161 PSF+Bulge+Disk 27.14 72.86 

NGC5247 PSF+Bulge+Disk 7.54 92.46 

NGC5371 PSF+Bulge+Disk+Bar 8.60 76.08 

NGC5427 PSF+Bulge+Disk 7.76 92.24 

NGC5483 Bulge+Disk+Bar 0.98 91.28 

NGC5643 PSF+Bulge+Disk+Bar 8.21 81.78 

NGC5676 PSF+Bulge+Disk 6.67 93.33 

NGC5701 PSF+Bulge+Disk+Bar 24.43 63.47 

NGC5713 PSF+Bulge+Disk 33.40 66.60 

NGC5850 PSF+Bulge+Disk+Bar 14.76 63.74 

NGC5921 PSF+Bulge+Disk+Bar 10.43 70.50 

NGC5962 Bulge+Disk+Bar 10.30 78.82 

NGC6215 PSF+Bulge+Disk 6.95 93.05 

NGC6221 Bulge+Disk+Bar 5.27 81.98 

NGC6300 PSF+Bulge+Disk+Bar 5.98 86.55 

NGC6384 Bulge+Disk 27.92 72.08 

NGC6753 PSF+Bulge+Disk+Bar 5.08 72.91 

NGC6782 PSF+Bulge+Disk+Bar 26.27 51.75 

NGC6902 PSF+Bulge+Disk 40.55 59.45 



14.20 



19.22 
11.35 
9.97 
14.98 



6.20 



19.06 



15.32 

7.73 
10.01 

12.09 

21.50 
19.06 
10.88 

12.76 
7.47 

22.02 
21.98 



0.66 
0.49 
0.11 
0.10 
0.45 
0.07 
0.39 
0.27 
0.18 
0.08 
0.11 
0.06 
0.68 
0.57 
0.37 
0.08 
0.11 
0.08 
0.01 
0.10 
0.07 
0.38 
0.50 
0.23 
0.15 
0.13 
0.07 
0.06 
0.07 
0.39 
0.07 
0.51 
0.68 



9.29 
21.43 
12.02 
2.82 
5.95 
4.31 
6.46 
9.52 
4.30 
4.14 
5.28 
5.53 
7.61 
4.68 
20.54 
8.46 
3.13 
4.74 
3.23 
5.68 
3.71 
11.13 
15.59 
6.22 
2.55 
2.33 
2.55 
3.16 
5.43 
15.66 
1.50 
3.90 
12.68 



0.67 
2.31 
1.04 
0.21 
0.51 
0.80 
1.14 
2.01 
0.25 
0.50 
0.65 
0.74 
0.96 
0.48 
3.38 
0.79 
0.55 
0.45 
0.39 
0.46 
0.54 
1.19 
2.02 
1.06 
0.26 
0.32 
0.27 
0.29 
0.42 
1.82 
0.32 
1.03 
2.43 



1.49 
1.75 
0.94 
0.47 
1.53 
0.85 
1.37 
2.78 
0.84 
0.05 
2.42 
1.48 
2.28 
2.41 
1.47 
1.04 
1.29 
0.60 
1.22 
2.14 
1.29 
2.41 
1.84 
1.58 
1.97 
1.37 
0.96 
2.73 
1.44 
3.13 
0.94 
1.10 
3.03 



58.16 
18.48 
32.82 
36.31 
30.61 
29.24 
45.02 
32.14 
23.23 
19.56 
52.65 
33.88 
28.39 
15.62 
43.98 
50.85 
55.81 
24.59 
23.25 
45.30 
23.10 
70.67 
18.52 
72.24 
40.88 
14.56 
14.02 
51.47 
47.08 
40.05 
19.17 
27.14 
30.13 



4.17 
1.99 
2.85 
2.69 
2.64 
5.45 
7.93 
6.79 
1.36 
2.34 
6.48 
4.55 
3.59 
1.61 
7.24 
4.78 
9.84 
2.31 
2.83 
3.63 
3.39 
7.56 
2.40 
12.29 
4.10 
1.99 
1.47 
4.78 
3.60 
4.65 
4.12 
7.20 
5.77 



52.58 



16.38 
17.07 
14.13 
32.96 



18.03 



37.85 



23.03 

8.74 
48.58 

26.01 

46.05 
35.05 
13.85 

19.33 
31.64 

7.88 
19.27 



3.77 



1.21 
1.47 
2.64 
5.80 



2.16 



4.78 



4.06 

1.06 
3.89 

2.78 

7.83 
3.51 
1.89 

1.80 
2.42 

1.69 
5.11 



0.50 



0.51 
0.44 
0.39 
0.38 



0.34 



0.31 



1.05 

0.32 
0.43 

0.40 

0.90 
0.92 
0.02 

0.87 
0.30 

0.80 
0.38 



o 











Table 2- 


—Continued 














Galaxy Name 


Best Fit 


B/T 


D/T 


Bar/T 


B/D 


Bulge Te 


Bulge re 


Bulge n 


Disk h 


Diskh 


Bar re 


Bar re 


Bar n 






% 


% 


% 




tl 


kpc 




// 


kpc 


// 


kpc 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


NGC6907 


PSF+Bulge+Disk+Bar 


10.32 


63.12 


26.56 


0.16 


3.49 


0.75 


1.10 


28.99 


6.19 


23.37 


4.99 


0.33 


NGC7083 


PSF+Bulge+Disk 


13.23 


86.77 


- 


0.15 


5.87 


1.23 


1.37 


23.26 


4.86 


- 


- 


- 


NGC7205 


PSF+Bulge+Disk 


6.08 


93.92 


- 


0.06 


4.70 


0.48 


0.88 


31.03 


3.16 


- 


- 


- 


NGC7213 


PSF+Bulge+Disk 


65.74 


34.26 


- 


1.92 


19.55 


2.42 


2.68 


63.19 


7.83 


- 


- 


- 


NGC7217 


PSF+Bulge+Disk 


53.71 


46.29 


- 


1.16 


21.07 


1.36 


2.21 


26.84 


1.73 


- 


- 


- 


NGC7412 


PSF+Bulge+Disk+Bar 


4.09 


68.60 


27.31 


0.06 


2.84 


0.34 


1.28 


23.58 


2.78 


35.06 


4.13 


1.20 


NGC7479 


PSF+Bulge+Disk+Bar 


8.76 


63.51 


27.72 


0.14 


6.00 


0.98 


1.09 


37.59 


6.17 


40.63 


6.67 


0.47 


NGC7552 


Bulge+Disk+Bar 


23.44 


61.08 


15.48 


0.38 


2.70 


0.31 


0.64 


17.54 


1.99 


42.05 


4.78 


0.24 


NGC7723 


PSF+Bulge+Disk+Bar 


5.11 


85.03 


9.85 


0.06 


2.61 


0.34 


0.54 


21.70 


2.82 


21.81 


2.84 


0.90 


NGC7727 


PSF+Bulge+Disk 


41.77 


58.23 


- 


0.72 


8.31 


1.05 


2.02 


24.46 


3.09 


- 


- 


- 


NGC7741 


Bulge+Disk+Bar 


3.09 


89.02 


7.90 


0.03 


9.51 


0.50 


0.31 


60.98 


3.18 


30.02 


1.56 


0.40 



Note. — Columns are: (1) Galaxy name. (2) The best fit chosen based on the criteria outlined in ij 13.31 (3) Bulge-to-total light ratio. If a PSF is included in 
the model, the contribution from the PSF is also added in. (4) Disk-to-total light ratio. (5) Bar-to-total light ratio. (6) Bulge-to-disk ratio (7) Bulge effective 
radius in arcseconds. (8) Bulge effective radius in kpc. (9) Bulge Sersic index. (10) Disk scalelength in arcseconds. (11) Disk scalelcngth in kpc. (12) Bar effective 
radius in arcseconds. (13) Bar Bulge effective radius in kpc. (14) Bar Sersic index. 
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Table 3. Decomposition For NGC 4643 



Fit Te or h Te or h n h/a Position Angle Fractional light 

(") (kpc) 

(1) (2) (3) (4) (5) (6) (7) (8) 



Stage 1 


Sersic 


27.90 


2.66 


4.44 


0.80 


-51.03 


100% 


Stage 2 


Bulge 


23.86 


2.30 


4.16 


0.80 


-51.08 


34.6% 




Disk 


335.88 


32.33 


1.00 


0.84 


66.94 


65.4% 


Stage 3 


Bulge 


5.43 


0.52 


2.53 


0.90 


60.52 


25.0% 




Disk 


48.22 


4.64 


1.00 


0.84 


66.94 


54.1% 




Bar 


21.30 


2.05 


0.62 


0.37 


-45.84 


20.9% 



Note. — Columns are: (1) Indicates whether the model is for Stage 1 (pure Sersic 
profile), Stage 2 (bulge+disk or bar+disk), or Stage 3 (bulge+disk+bar). (2) One 
of pure Sersic profile, Sersic bulge, exponential disk, or Sersic bar. (3) Effective 
radius of the bulge/bar or disk scalelength in arcseconds. (4) Effective radius of the 
bulge/bar or disk scalelength in kpc. (5) The Sersic index of the profile. (6) Axis 
ratio of the component. (7) Position angle in degrees of the component. Positive 
angles are East of North. (8) The percentage of total galaxy light contributed by 
the component. 
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Table 4. Decomposition For NGC 4548 



Fit Te or h Te or h n h/a Position Angle Fractional light 

(") (kpc) 

(1) (2) (3) (4) (5) (6) (7) (8) 



Stage 1 


Sersic 


154.59 


5.19 


5.19 


0.80 


78.31 


100% 


Stage 2 


Bulge 


57.86 


1.94 


4.32 


0.76 


75.77 


61.5% 




Disk 


60.39 


2.03 


1.00 


0.75 


-32.54 


38.5% 


Stage 3 


Bulge 


6.98 


0.23 


1.56 


0.88 


-66.50 


13.0% 




Disk 


58.22 


1.96 


1.00 


0.75 


-32.54 


68.6% 




Bar 


44.91 


1.51 


0.51 


0.35 


66.65 


18.4% 



Note. — Columns are: (1) Indicates whether the model is for Stage 1 (pure Sersic 
profile), Stage 2 (bulge-|-disk or bar-|-disk), or Stage 3 (bulge-|-disk-|-bar). (2) One 
of pure Sersic profile, Sersic bulge, exponential disk, or Sersic bar. (3) Effective 
radius of the bulge/bar or disk scalelength in arcseconds. (4) Effective radius of the 
bulge/bar or disk scalelength in kpc. (5) The Sersic index of the profile. (6) Axis 
ratio of the component. (7) Position angle in degrees of the component. Positive 
angles are East of North. (8) The percentage of total galaxy light contributed by 
the component. 
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Table 5. Decomposition For NGC 4902 



Fit Te or h Te or h n h/a Position Angle Fractional light 

(") (kpc) 

(1) (2) (3) (4) (5) (6) (7) (8) 



Stage 1 


Sersic 


154.4 


28.8 


5.11 


0.54 


69.0 


100% 


Stage 2 


Bulge 


11.5 


2.15 


1.02 


0.38 


68.5 


17.4% 




Disk 


32.8 


6.12 


1.00 


0.84 


81.1 


82.6% 


Stage 3 


Bulge 


4.31 


0.80 


0.85 


0.68 


-52.96 


5.59% 




Disk 


29.2 


5.45 


1.00 


0.84 


81.1 


82.6% 




Bar 


14.1 


2.64 


0.39 


0.22 


66.37 


9.97% 



Note. — Columns are: (1) Indicates whether the model is for Stage 1 (pure Sersic 
profile), Stage 2 (bulge+disk or bar+disk), or Stage 3 (bulge+disk+bar). (2) One 
of pure Sersic profile, Sersic bulge, exponential disk, or Sersic bar. (3) Effective 
radius of the bulge/bar or disk scalelength in arcseconds. (4) Effective radius of the 
bulge/bar or disk scalelength in kpc. (5) The Sersic index of the profile. (6) Axis 
ratio of the component. (7) Position angle in degrees of the component. Positive 
angles are East of North. (8) The percentage of total galaxy light contributed by 
the component. 
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Table 6. Checking GALFIT robustness with different input guesses for Stage 3 

B/T Bulge re Bulge n D/T Disk h Bar/T Bar r^ Bar n 



Disk h 


Bar/T 


Bar r, 


(") 




{") 


(6) 


(7) 


(8) 



_« (2) (3) (4) (5) (6) (7) (8) (9) 

NGC 4548 



Initial guesses from ID decomposition 


17.5% 


7.39 


1.17 


63.5% 


28.4 


19.0% 


37.5 


0.54 


Stage 3 Output 


13.0% 


6.98 


1.56 


68.6% 


58.2 


18.4% 


44.9 


0.51 


Initial guesses from from Stage 2 


11.1% 


7.50 


1.70 


69.9% 


64.5 


19.1% 


37.5 


0.54 


Stage 3 Output 


13.0% 


6.98 


1.56 


68.6% 


58.2 


18.4% 


44.9 


0.51 






NGC 


4643 












Initial guesses from ID decomposition 


33.6% 


7.18 


0.86 


40.4% 


37.5 


26.0% 


22.0 


0.60 


Stage 3 Output 


25.0% 


5.43 


2.53 


54.1% 


48.2 


20.9% 


21.3 


0.62 


Initial guesses from Stage 2 


24.1% 


5.30 


2.5 


51.8% 


46.4 


24.1% 


22.0 


0.60 


Stage 3 Output 


25.0% 


5.43 


2.53 


54.1% 


48.2 


20.9% 


21.3 


0.62 



Note. — For each galaxy, the rows are: (1) The initial guesses to Stage 3 provided by ID decomposition. (2) The 
corresponding output for the input of row (1). (3) The initial guesses to Stage 3 provided by Stage 2. (4) The corresponding 
output for the input of row (3). 

The columns are: (2) Bulge-to-total luminosity fraction. (3) Effective radius of the bulge in arcseconds. (4) The Sersic 
index of the profile. (5) Disk-to-total luminosity fraction. (6) Disk scalelength in arcseconds. (7) Bar-to-total luminosity 
fraction. (8) Bar effective radius in arcseconds. (9) Bar Sersic index. 
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Table 7. Parameter Coupling 



Type of Fit 


X^ 


Bulge Te 

(") 


Bulge n 


B/T 


Disk h 


D/T 


Bar re 


Bar n 


Bar/T 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 










NGC 


3885 












initial input 


11241.33 


2.32 


2.00 


27.7% 


10.44 


72.3% 


- 


- 


- 


n=free 


8210.12 


3.13 


0.46 


27.1% 


11.94 


72.9% 


- 


- 


- 


n=l fixed 


[8307.82] 


[2.96] 


[1.00] 


[29.1%] 


[12.90] 


[70.9%] 


- 


- 


- 


n=4 fixed 


[7373.55] 


[6.28] 


[4.00] 


[53.7%] 


[20.45] 


[46.3%] 


- 


- 


- 


n=l re-fit 


8210.12 


3.13 


0.46 


27.1% 


11.94 


72.9% 


- 


- 


- 


n=4 re-fit 


7258.46 


7.11 


4.49 


57.3% 


19.87 


42.7% 


- 


- 


- 








NGC 


7213 












initial input 


605757.39 


5.80 


4.00 


56.4% 


9.25 


43.6% 


- 


- 


- 


n=free 


71170.10 


19.55 


2.68 


65.7% 


63.19 


34.3% 


- 


- 


- 


n=l fixed 


[96312.29] 


[7.81] 


[1.00] 


[31.4%] 


[27.79] 


[68.6%] 


- 


- 


- 


n=4 fixed 


[70998.69] 


[41.11] 


[4.00] 


[95.4%] 


[6.75] 


[4.6%] 


- 


- 


- 


n=l re-fit 


71170.10 


19.53 


2.68 


65.7% 


63.16 


34.3% 


- 


- 


- 


n=4 re-fit 


69970.51 


72.04 


5.60 


91.9% 


7.98 


8.1% 


- 


- 


- 








NGC 


4643 












initial input 


7841.46 


5.34 


2.50 


24.1% 


46.40 


51.8% 


22.04 


0.60 


24.1% 


n=free 


2111.59 


5.43 


2.53 


25.0% 


48.22 


54.1% 


21.30 


0.62 


20.9% 


n=l fixed 


[3001.98] 


[2.31] 


[1.00] 


[9.6%] 


[53.37] 


[57.0%] 


[15.05] 


[0.88] 


[33.4%] 


n=4 fixed 


[2491.06] 


[12.08] 


[4.00] 


[35.4%] 


[87.22] 


[51.7%] 


[23.37] 


[0.61] 


[12.9%] 


n=l re- fit 


2111.59 


5.43 


2.53 


25.0% 


48.21 


54.1% 


21.30 


0.62 


20.9% 


n=4 re-fit 


2111.59 


5.43 


2.53 


25.0% 


48.22 


54.1% 


21.30 


0.62 


20.9% 


NGC 4151 


initial input 


1388891.11 


7.50 


0.50 


24.6% 


33.00 


65.0% 


55.50 


0.10 


10.3% 


n=free 


16391.66 


9.36 


0.44 


41.3% 


37.74 


50.3% 


57.00 


0.10 


8.4% 


n=l fixed 


[18100.67] 


[11.58] 


[1.00] 


[42.2%] 


[86.28] 


[41.1%] 


[56.70] 


[0.17] 


[16.7%] 


n=4 fixed 


[20051.12] 


[108.39] 


[4.00] 


[51.4%] 


[8.21] 


[20.4%] 


[68.23] 


[0.28] 


[28.2%] 


n=l re- fit 


16391.58 


9.38 


0.44 


41.2% 


38.31 


49.7% 


55.57 


0.12 


9.0% 


n=4 re-fit 


16390.93 


9.36 


0.44 


41.3% 


38.20 


49.8% 


55.88 


0.11 


8.8% 



Note. — For each galaxy, the rows are: (1) The model outputs (columns 2 to 10) from GALFIT based on the 
Stage 1 or Stage 2 fits. These outputs are used as initial guesses in the fits for rows 2 and 3. (2) The model 
outputs (columns 2 to 10) from GALFIT when the bulge n is allowed to vary freely, and the input guesses are 
based on the parameters in columns 2 to 10 of row 1. (3) The model outputs (columns 2 to 10) from GALFIT 
when the bulge Sersic index is held fixed at n = 1, and the input guesses are based on the parameters in columns 
2 to 10 of row 1. The model is not always physically meaningful because the bulge index is fixed. (4) As in row 
3, except that the bulge Sersic index is now held fixed at n = 4. The model is not always physically meaningful 
because the bulge index is fixed. (5) The model outputs (columns 2 to 10) from GALFIT, when the bulge n is 
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allowed to vary freely, but the input guesses are now based on the parameters in columns 2 to 10 of row 3. (6) 
The model outputs (columns 2 to 10) from GALFIT, when the bulge n is allowed to vary freely, but the input 
guesses are now based on the parameters in columns 2 to 10 of row 4. The columns are: (1) The type of fit whose 
model outputs are shown in columns 2 to 10. (2) The x^ of the fit. (3) Effective radius of the bulge in arcseconds. 
(4) The Sersic index of the bulge (in the case of Stage 2 bulge-disk fits in row 1 and in the case of all fits in 
rows 2 to 6), or the Sersic index of a single component (in the case of Stage 1 fits in row 1). (5) Bulge-to-total 
luminosity fraction. (6) Disk scalelength in arcseconds. (7) Disk-to-total luminosity fraction. (8) Bar effective 
radius in arcseconds. (9) Bar Sersic index. (10) Bar-to-total luminosity fraction. 
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Table 8. Mass Breakdown of Galactic Structures 



Structure Mass 



Bulges 


18.9% ± 3.49% 


Disks 


71.6% ±4.02% 


Bars 


9.58% ± 2.62% 



Bulges with n > 2 10.4% ± 2.72% 
Bulges with n < 2 8.42% ± 2.47% 
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Table 9. B /T and bulge n in bright/high mass spirals 

Sample {Mb < -19.3) 

of spirals with bulge n > 4 1.40% ± 0.98% 

of spirals with bulge 2 < n < 4 22.4% ± 3.49% 

of spirals with bulge n < 2 76.2% ± 3.56% 

of spirals with bulge n>2 23.8% ± 3.56% 

of spirals with B/T < 0.2 68.5% ± 3.88% 

of spirals with B/T > 0.2 31.1% ± 3.88% 

of spirals with 0.2 < B/T < 0.4 23.1% ± 3.52% 

of spirals with B/T > 0.4 8.39% ± 2.32% 



Sample M^ > 1 x lO^^M, 







of spirals with bulge n > 4 1.77% ± 1.24% 

of spirals with bulge 2 < n < 4 23.9% ± 4.01% 

of spirals with bulge n < 2 74.3% ±4.11% 

of spirals with bulge n>2 25.7% ±4.11% 

of spirals with B/T < 0.2 66.4% ± 4.44% 

of spirals with B/T > 0.2 33.6% ± 4.44% 

of spirals with 0.2 < B/T < 0.4 25.7% ± 4.11% 

of spirals with B/T > 0.4 7.96% ± 2.55% 
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Table 10. Bar Fraction as a Function of B /T and Bulge Index 



Sample {Mb < -19.3) 



ion 
ion 
ion 
ion 



Bar fract 
Bar fract 
Bar fract 
Bar fract 
Bar fraction 
Bar fraction 
Bar fraction 
Bar fraction 



in spirals with bulge n > 4 
in spirals with bulge 2 < n < 4 
in spirals with bulge n < 2 
of spirals with bulge n > 2 
of spirals with B/T < 0.2 
of spirals with B/T > 0.2 
of spirals with 0.2 < B/T < 0.4 
of spirals with B/T > 0.4 



0.00% ± 0.00% 
37.5% ± 8.56% 
65.1% ±4.56% 
35.3% ± 8.20% 
68.4% ±4.70% 
35.6% ±7.14% 
42.4% ± 8.60% 
16.7% ± 10.8% 



Sample M^ > 1 x lO^^M, 







Bar fraction of spirals with bulge n > 4 
Bar fraction of spirals with bulge 2 < n < 4 
Bar fraction of spirals with bulge n <2 
Bar fraction of spirals with bulge n > 2 
Bar fraction of spirals with B/T < 0.2 
Bar fraction of spirals with B/T > 0.2 
Bar fraction of spirals with 0.2 < B/T < 0.4 
Bar fraction of spirals with B/T > 0.4 



0.00% ± 0.00% 
44.4% ± 9.56% 
63.1% ±5.27% 
41.4% ±9.15% 
68.0% ± 5.39% 
36.8% ± 7.83% 
41.4% ±9.15% 
22.2% ± 13.9% 
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Table 11. B /T: Data versus Hierarchical Models of Galaxy Evolution 



(1) 



Data Model Model spirals with Model spirals with All model 

major mergers no major merger spirals 
since redshift z\ since redshfit z\ 

(2) (3) (4) (5) (6) 



Fraction of spirals with B/T < 0.2 66.4% ± 4.44% SK 

Fraction of spirals with 0.2 < B/T < 0.4 25.7% ± 4.11% SK 

Fraction of spirals with 0.4 < B/T < 0.75 8.0% ± 2.5% SK 

Fraction of spirals with B/T < 0.75 100% SK 



For zi <2 
0.1% 
2.2% 
11.2% 
13.5% 



For zi <2 
67.2% 
16.6% 

2.6% 
86.5% 



67.3% 
18.8% 
13.8% 
100.0% 



Fraction of spirals with B/T < 0.2 66.4% ± 4.44% SK 

Fraction of spirals with 0.2 < B/T < 0.4 25.7% ± 4.11% SK 

Fraction of spirals with 0.4 < B/T < 0.75 8.0% ± 2.5% SK 

Fraction of spirals with B/T < 0.75 100% SK 



For zi <4 
1.6% 
5.5% 
12.6% 
19.8% 



For 21 < 4 
65.7% 
13.3% 
1.3% 
80.2% 



67.3% 
18.8% 
13.9% 
100.0% 



Fraction of spirals with B/T < 0.2 66.4% ± 4.44% PH 

Fraction of spirals with 0.2 < B/T < 0.4 25.7% ± 4.11% PH 

Fraction of spirals with 0.4 < B/T < 0.75 8.0% ± 2.5% PH 

Fraction of spirals with B/T < 0.75 100% PH 



For zi <2 
3.0% 
7.1% 
5.9% 
15.7% 



For 21 < 2 
67.7% 
13.1% 
3.1% 
84.0% 



70.7% 
20.1% 
9.0% 



Note. — Column 2 shows the empirical fraction of high mass spirals with a given B/T. Columns 3 to 6 show the fraction of 
high mass spirals with a given B/T in the models, which define spirals as systems with B/T < 0.75. Column 3 indicates the model 
used: SK refers to the models by Khochfar & Burkert (2005) and Khochfar & Silk (2006), while PH refers to those by Hopkins et 
al. (2009). The models are described in ^5.8l and have different prescriptions for the amount of mass that gets converted into bulge 
stars. Column 4 shows the fraction of spirals, which have a given B/T and have had one or more major mergers since redshift zi. 
Column 5 shows the fraction of spirals, which have a given B/T and have had no major merger since redshift zi . Column 6 is the 
sum of column 4 and 5, and shows the fraction of spirals with a given B/T in the models. It can be directly compared to the data 
in column 2. 



-62- 



Table 12. Comparison of model merger fraction since z < 4 and bulge Sersic index 



Stellar mass limit: 1.0 xlO^^Mp 







Model galaxies: 7674 
Observed galaxies: 113 



Test 1 



% of model galaxies w/ B/T < 0.2 with no major merger since z < 4: 5040/7674 = 65.7% ± 0.54% 

% of observed galaxies w/ B/T < 0.2 and bulge n < 3: 75/113 = 66.4% ± 4.44% 

% of observed galaxies w/ B/T < 0.2 and bulge n < 2.5: 74/113 = 65.5% ± 4.47% 

% of observed galaxies w/ B/T < 0.2 and bulge n < 2: 67/113 = 59.3% ± 4.62% 



Test 2 



% of model galaxies w/ B/T > 0.4 and a major merger since z < 4: 
% of observed spirals w/ B/T > 0.4 and bulge n > 3: 
% of observed spirals w/ 5/T > 0.4 and bulge n > 2.5: 
% of observed spirals w/ 5/T > 0.4 and bulge n > 2: 



971/7674 = 12.7% ± 0.38% 
2/113 = 1.77% ± 1.24% 
4/113 = 3.54% ±1.74% 
7/113 = 6.19% ± 2.27% 



Test 3 



% of model galaxies w/ B/T < 0.2 and a major merger since z < 4: 
% of observed galaxies w/ B/T < 0.2 and bulge n > 3: 
% of observed galaxies w/ B/T < 0.2 and bulge n > 2.5 
% of observed galaxies w/ i?/T < 0.2 and bulge n > 2 



126/7674 = 1.64% ±0.15% 
0/113 = 0.00% ±0.00% 
1/113 = 0.88% ±0.88% 
8/113 = 7.08% ±2.41% 
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Fig. 1. — Our final sample SI consists of 143 bright {Mb < —19.3) low-to-nioderately 
inclined {i < 70°) spirals in the OSUBSGS survey. The distribution of absolute 5-band 
magnitude for the sample of bright spirals in the OSUBSGS survey is shown in the top panel 
before (solid line) and after (shaded greyscale) the cut to remove highly inclined (i > 70°) 
spirals. The distribution of Hubble types for the sample is shown in the bottom panel before 
(solid line) and after (shaded greyscale) the cut to remove highly inclined {i > 70°) spirals. 
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Fig. 2. — The luminosity function of the full OSUBSGS sample is compared with the 5-band 
Schechter luminosity function (SLF). The former is calculated as described in § 12.11 using 
equation ([1]). The parameters for the SLF are $* = 5.488 x 10"'^ Mpc~'^, a = —1.07, and 
M^ = -20.5 (Efstathiou, EUis & Peterson 1988), corresponding to Ho=70 km/s Mpc'^ 
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Fig. 3. — Out of our final sample SI of 143 bright {Mb < —19.3) low-to-moderately 
inclined (z < 70°) OSUBSGS spirals, stellar masses could be estimated for 126 galaxies. 
Their stellar mass distribution is shown, as determined in § 12. 2[ Most have stellar masses 
M^ > 1.0 X lO^'' Mq. This sample of 126 galaxies is referenced henceforth as the sample S2. 
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Stage 1: 
One Component 



Input guesses for 
single Sersic component 



Fit single Sersic profile 



Stage 2: 

Two Components 

(bulge + disk or bar+disk) 



Use Stage 1 output as 

guesses for input bulge parameters 

(disk b/a, PA fixed to pre-determined values) 



Fit Sersic + exponential profile 



Stage 3: 

Three Components 

(bulge+disk+bar) 



Use Stage 2 output as 

guesses for bulge and disk parameters; 

input guesses for bar parameters 

(disk b/a, PA fixed to pre-determined values) 



Fit Sersic bulge 

+ exponential disk 

+ Sersic bar 



Choose the best fit from Stage 2 and Stage 3 based on: 

Residual map 

Model parameters 

Data image 



Fig. 4. — An overview of the method of decomposition. All images are subjected to Stages 
1, 2, and 3. Either the best fit of Stage 2 or Stage 3 is chosen as the best model. 
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Bulge + disk + bar 



Fig. 5. — Complete ID decomposition for NGC 4643. Note the prominent bar residuals in 
the residual for the Stage 1 and Stage 2 bulge-disk decomposition. This is a case where the 
prominent bar causes the Stage 2 bulge-disk fit to artificially extend the bulge and inflate 
the BjT . The disk fitted in Stage 2 has a low surface brightness and is very extended, 
well beyond the real disk: the hja and PA of the fitted disk is shown as contours. Stage 3 
bulge-disk-bar decomposition provides the best model. The x^ for the Stage 1, Stage 2, and 
Stage 3 residual images are 7360.7, 7284.8, and 2111.59, respectively. See Table [3] for the fit 
parameters. 
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Fig. 6. — The complete 2D decomposition for NGC 4548. This is an extreme example where 
the prominent bar results in an extended bulge and inflated B /T in the Stage 2 bulge-disk 
fit. Like NGC 4643 in Figure O the disk fitted in Stage 2 has a low surface brightness and is 
very extended: its h/a and PA are shown as contours. Stage 3 bulge-disk-bar decomposition 
provides the best model. The Xv for the Stage 1, Stage 2, and Stage 3 residual images are 
7076.1, 6301.3, and 3260.4, respectively. See Table H] for the fit parameters. 
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Data 




5 kpc H-band 20- 

Fig. 7. — This plot shows the data image, Stage 2 model, and Stage 3 model for NGC 4902. 
The Stage 2 bulge is too bright and is extended along the major axis of the bar {B/T=31.2% 
and b/a=0A5). In Stage 3, the bulge and bar are fit with distinct components {B/T=5.59%, 
bulge 6/a=0.68, Bar/T=9.97%, bar 6/a=0.22). All other fit parameters appear in Table O 
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Fig. 8. — The data images and Stage 3 bulge-disk-bar decomposition models of NGC 5427 
and NGC 7412 are shown. The Stage 3 models each distinctly show a false bar component, 
which is not present in the data images. The false components can be inspired by prominent 
spiral arms, such as those present in these galaxies. Such cases are flagged during the visual 
inspection of fits and the Stage 3 bulge-disk-bar decomposition is discarded in favor of the 
Stage 2 bulge-disk decomposition. 
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Fig. 9. — This plot compares the B/T obtained by fitting the bulge of NGC 4548 with a 
Sersic model of constant h/a as opposed to a Sersic model varying h/a. To mimic a Sersic 
model with varying h/a in GALFIT, the bulge was fitted with ten concentric Sersic profiles 
with fixed re, each separated by 0.75". The top two panels show the run of h/a and PA of 
the ten concentric Sersic profiles. The bottom panel shows the cumulative B/T calculated 
by summing all models with r < r^.. The bulge h/a (0.88), PA (-66.5), and B/T (13%) from 
the original Sersic fit of constant h/a (Table H]) are indicated with horizontal lines on the 3 
panels. 
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Fig. 10. — An elementary test is to determine if GALFIT can recover the known parameters 
of artificial noisy images. Noisy images were simulated by taking parametric model images 
(left panels) produced by GALFIT, and adding noise and sky background (right panels). 
The noisy images were then fitted to see if the original known parameters can be recovered. 
See § 14.21 for details. 
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Fig. 11. — The plots compare recovered versus original model parameters for the simulated 
images discussd in § 14.21 The vertical axis limits demonstrate the range explored for each 
parameter. The dotted line shows y = x for comparison. Except for some extreme cases 
where the images were highly distorted by noise, all parameters were recovered to within a 
few percent. 
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Fig. 12. — The ratio of model-to-recovered parameter is plotted against mean surface 
brightness inside the disk scalelength, yU = mag + 2.5 x Zof7io(2 x vr x h/a x /i^), for the 
simulated images discussed in § 14.21 Surface brightness is not photometrically calibrated 
and is shown for a zeropoint of 0. 
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Fig. 13. — The top, middle, and bottom panels show stellar mass for bulges, disks, and bars, 
respectively, along the Hubble sequence. Values are shown for sample S2 of 126 galaxies in 
Fig. [31 The legend in each panel indicates the type of decomposition used for each data 
point. 
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Fig. 14. — The individual and mean B/T (left panels) and bulge Sersic index (right panels) 
are plotted as a function of Hubble type for the sample SI of bright galaxies, and as a function 
of galaxy stellar mass for sample S2. The error bars indicate the standard deviation of the 
population around the mean in each bin. The legend in each panel indicates the type of 
decomposition used for each data point. The mean B/T and bulge index in barred galaxies 
differ systematically from unbarred galaxies, but there is a large overlap in the individual 
values. As many as ~ 69% of bright spiral galaxies have B/T < 0.2; these bulges are 
pervasive and exist across the Hubble sequence. Furthermore, as many as ~ 76% of bright 
spirals have low n < 2 bulges. Such bulges exist in barred and unbarred galaxies across a 
wide range of Hubble types. 
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Fig. 15. — The relation between B /T and bulge index is shown. In the top panel, galaxies 
are coded according to bar class. The legend indicates the type of decomposition used for 
each data point. In the lower panel, galaxies are coded according to Hubble type. A striking 
~ 76% of bright spirals have low n < 2 bulges. Such bulges exist in barred and unbarred 
galaxies across a wide range of Hubble types, and their B /T range from 0.01 to 0.4, with 
most having B/T < 0.2. A moderate fraction (~ 22%) have intermediate 2 < n < 4 bulges. 
These exist in barred and unbarred SO/a to Sd galaxies, and their B/T spans a wide range 
(0.05 to 0.5). Only (~ 1%) have n > 4. 
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Fig. 16. — B/T is plotted against Bar/T and sorted by bulge Sersic index. There are six 
galaxies with Bar/T > 0.3. 
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Fig. 17. — The top two rows show //-band images of barred galaxies, which have early 
RC3 Hubble types, but yet have BjT < 0.2. The bottom row shows iJ-band images of 
unbarred galaxies, which have late RC3 Hubble types, but yet have BjT ~ 0.4. The Hubble 
types assigned to these objects more reflect disk smoothness and spiral arm topology than 
BjT . All images are from OSUBSGS with characteristics as described in § [21 
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Fig. 18. — B/D is plotted against Hubble type. The legend indicates the type of decompo- 
sition used for each data point. The mean values for barred and unbarred together in each 
bin are shown. 
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Fig. 19. — The properties of bars are shown. The error bars indicate the standard deviation 
of the population around the mean in each bin. The legend in each panel indicates the type 
of decomposition used for each data point. Upper left: Mean and individual Bar/T plotted 
against Hubble type. Upper right: Mean and individual bar Sersic indices plotted against 
Hubble type. Lower left: Bar/T plotted against total galaxy stellar mass. The mean Bar/T 
in bins of stellar mass is indicated. Lower right: Bar Sersic index plotted against total galaxy 
stellar mass. 
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Fig. 20. — Bar strength is plotted against Hubble type and the properties of bars are 
compared with bulges. The legend in each panel indicates the type of decomposition used 
for each data point. Upper left: Bar strength, the product of Bar/T and peak bar ellipticity 
Char from MJ07 is plotted against Hubble type. Upper right: Bar/T is plotted against bar 
Sersic index. Lower left: Bar/T is plotted against bulge Sersic index. Lower right: Bar/T is 
plotted against B/T. In the first plot, mean bar strong is calculated for each Hubble type. 
In the latter three plots, mean Bar/T is calculated for bins along the ordinate axis. The 
error bars indicate the standard deviation of the population around the mean in each bin. 
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Fig. 21. — iJ-band OSUBSGS images of spirals with prominent bars of Bar/T (> 0.3) are 
shown. An interesting example is the oval or lens galaxy NGC 1317: the bar has a low 
ellipticity, but its B /T is large because it is extended and massive. Such bars/lenses may 
exert significant gravitational torques although they are not very elongated. 
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Fig. 22. — Shown is the distribution of Sersic indices n for remnants of 1:1 gas-rich major 
mergers in the simulations of Hopkins et ah 2008: they he in the range of 2 < n < 6. 
Specifically, ~ 22% of the remnants have classical n > 4, as much as 20% have low n < 2.5, 
while 50% have n < 3. Almost none have n < 2. [Figure: courtesy of Phil Hopkins] 
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Fig. 23. — For those high mass (M^ > 1.0 x lO^^M©) galaxies in the theoretical models that 
experienced a major merger (see § 15.81) . the B/T of the remnant at z ~ is plotted against 
the redshift 2;iast of the last major merger. Systems where the last major merger occurred at 
earlier times have had more time to grow a disk and thus have a lower B/T at 2 ~ 0. The 
dispersion in the present-day B/T at a given 2;iast is due to the different times spent by a 
galaxy in terms of being a satellite versus a central galaxy in a DM halo, since the cooling of 
gas and the growth of a disk is stopped when a galaxy becomes a satellite. In the model, a 
high mass galaxy that has undergone a major merger at z < 2 has a present-day B/T > 0.2. 
In effect, a high mass spiral can have a present-day B/T < 0.2 only if its last major merger 
occurred at 2; > 2 (lookback times > 10 Gyr). 
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Fig. 24. — The cumulative fraction F of high mass (M^ > 1.0 x 10^'' Mq) spirals with 
present-day B /T < the x-axis value is shown for the data (colored lines/ points) and for the 
theoretical model (black lines/points) described in § 15.81 Model and data spirals are defined 
as systems with B/T < 0.75. The magenta line shows F from the data, while the other two 
colored lines break this F in terms of bar class (top panel) or bulge n (lower panel). The 
black dashed line shows F from all model galaxies, while the black dotted line and black 
dots show the contribution of model galaxies that experienced, respectively, no major m,erger 
and one or more major mergers since z < 4. Major mergers are defined here as those with 
M1/M2 > 1/4. In the model, the fraction (~ 1.6%; see Table fTTl) of high mass spirals, which 
have undergone a major merger since z < 4 and host a bulge with a present-day B/T < 0.2 
is a factor of over 20 smaller than the observed fraction (~ 66%) of high mass spirals with 
present-day B/T < 0.2. Thus, bulges built via major mergers since z < 4 seriously fail to 
account for most of the low B/T < 0.2 bulges present in two-thirds of high-mass spirals. 
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Fig. 25. — This figure is similar to Fig. [24l except that the model now defines major mergers 
as those with mass ratio M1/M2 > 1/6. 
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Fig. 26. — This figure is similar to Fig. [211 except that here spirals are considered to be 
systems with a B/T < 0.55 rather than 0.75 in the models, and a corresponding cut is 
applied to the data points. 
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Fig. 27. — This figure is similar to Fig. [2U except that BjT of all the observed galaxies has 
been multiplied by a factor of two, in order to test what would happen in the case where the 
MjL ratio of the bulge in if -band is twice as high as that of the disk and bar. This could 
happen in an extreme example where the dominant bulge stellar population was much older 
(e.g. 12 Gyr) than the age of the dominant disk stellar population (e.g., 3 Gyr). In such a 
case, the fraction of high mass spirals with B jT < 0.2 would change from ~ 66% in Fig. [2l] 
to ~ 50%. However, the main conclusion that bulges built by major mergers since z < A 
cannot account for most of the low B/T < 0.2 bulges, present in a large percentage (~ 55%) 
of spirals still holds. 



